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I .  INTRODUCTION  ^ 

A.  Purpose 

A program  has  been  written  to  incorporate  the  necessary  routines 
required  to  provide  more  expedient  computer  plotting  capability  for 
general  purpose  two  dimensional  and  three  dimensional  (x,y)  data 
plotting.  The  intent  of  this  program,  INTERPOL*,  is  to  minimize 

programming  effort.  The  interactive  format  of  the  program  reduces  the  ! 

typical  three-step  planning,  writing,  and  debugging  process  to  an  i 

efficient,  single-step  question  and  answer  execution.  Foreknowledge  of  j 

computer  plotting  techniques  is  not  necessary,  although  some  experience  I 

may  be  helpful. 

B.  Scope 

INTERPOL  was  written  for  use  on  the  UNIVAC  1108  computer.  The  plot- 
ting routines  employed  are  standard  Cal  Compl-3  FORTRAN  subprograms. 

Driver  routines  and  miscellaneous  I/O  routines  all  use  standard  FORTRAN. 

Conversion  of  INTERPOL  to  other  machines  employing  off-line  Cal  Comp 
pen  plotters  can  be  accomplished  with  little  effort. 

■ 

Use  of  this  program  is  most  advantageous  with  CRT  time-sharing 
terminals,  but  TTY  type  terminals  may  be  used.  Length  of  execution  is 
dependent  on  time -sharing  response,  number  of  graphs  generated,  and 
amount  of  data  processed  per  graph  (the  former  of  these  being  dominant) . 

INTERPOL  requires  -vSOK  of  UNIVAC  1108  core. 


II.  PROGRAM  STRUCTURE 

A.  Problem  Flow 

Figure  1 shows  the  basic  flow  structure  for  INTERPOL.  The  user 
is  first  asked  for  the  logical  unit  number  upon  which  the  digital 
graphic  commands  will  be  written.  The  pen  plotter  is  then  initialized 


*Interaative  jotting  for  Off-Line  Systems 

1,  'Programming  Cal  Comp  Pen  Plotters" ^ California  Computer  Products 
Phamphlet  lOOSAj  Anaheim^  California^  September  1969. 

2,  "User's  Manual",  Cal  Comp  Graphics  Functional  Software  - USAS 
FORTRAN  General,  California  Computer  Products,  Anaheim,  California, 
November  1968. 

3,  "User's  Manual"  Cal  Comp  Graphics  Functional  Software  -USAS 
FORTRAN /Scientific,  California  Computer  Products,  Anaheim 
California,  July  1969. 
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Figure  1.  INTERPOL  Flow  Diagram 


I 


accordingly  and  problem  execution  begins, 
loops;  one  controlling  the  general  design 
or  three-dimensional,  title,  axis  labeling 
other  controlling  the  display  of  the  data 
or  line/symbol  combination).  The  user  has 
any,  or  all,  input  parameters  within  each 
design  of  successive  graphs  or  representat 
In  this  manner  INTERPOL  allows  for  multipl 
and  multiple  graphs  during  one  execution. 

B.  Cal  Comp  Routines 


Each  run  is  comprised  of  two 
of  the  graph  (two  dimensional 
, scaling,  and  grid);  and  the 
(type  of  line,  type  of  symbol, 
the  opportunity  to  change 
loop,  thereby  changing  the 
ion  of  the  data  sets  at  will, 
e sets  of  data  on  one  graph 


Table  I lists  the  Cal  Comp  routines  employed  by  INTERPOL  along  with 
a brief  description  of  the  function  they  perform.  The  arguments  used 
in  the  call  to  these  routines  are  all  set  internally,  based  on  the  user's 
selection  of  the  graphic  design  parameters.  The  routines  provided  are 
among  the  most  generalized  available  from  Cal  Comp.  Specialized  sub- 
programs which  draw  specific  figures  or  curves  were  intentionally  omit- 
ted in  order  to  keep  INTERPOL  as  efficient  as  possible  for  plotting 
(x,y)  data  only.  Thus,  the  user  has  available  to  him  linear,  semi- 
logarithmic,  or  log-log  scaling,  with  smooth  line  fit,  straight  line 
point  connection,  dashed  line  point  connection,  symbols  printed  at 
points  only,  or  line/symbol  combinations.  This  is  adequate  for  most 
graph  applications  (examples  in  Appendix  A)  and  keeps  the  user-supplied 
input  at  a minimum. 

C.  INTERPOL  Routines 

Table  II  lists  the  routines  used  by  INTERPOL  to  provide  the  inter- 
face between  the  user  and  the  Cal  Comp  software.  Brief  descriptions 
are  given  in  the  table  but  several  require  further  explanation. 


TITLE 


NUGRID  -- 


- This  subroutine  prints  a title  block  if  requested  by  the  user. 
It  provides  up  to  five  lines  of  script,  25  characters  per  line. 
The  user  may  position  the  title  block  anywhere  on  the  graph  by 
specifying  an  upper  left  coordinate  (i.e.,  where  the  first  line 
is  to  begin)t  At  its  maximum,  the  title  block  size  is  approxi- 
mately .3.5  inches  horizontally  by  1.5  inches  vertically. 

--  This  subroutine  draws  either  a vertical,  horizontal  or  vertical 
and  horizontal  grid  pattern  on  the  plot.  The  increment  of  each 
line  along  an  axis  is  dictated  by  the  scaling  of  the  data  (tick 
marks)  for  that  axis.  NUGRID  is  compatible  with  TITLE  in  that 
the  lines  drawn  will  not  intersect  the  script  of  the  title 
block. 


*This  coordinate  is  an  (x,y)  (joint,  in  inches,  relative  to 
the  origin  of  the  graph. 


n 


Table  I.  Standard  Cal  Comp  Routines  Used  by  INTERPOL 


N^ME 

PIXITS 

PLOT 

FACTOR 

SCALE 

AXIS 

LC.AXIS 

SYMBOL 

LINE 

DASHLN , 

FLINE 

LGLIN 


FUNCTION 

Initializes  the  pen  plotter. 


Establishes  origin:  draws  a line. 


Expands  or  reduces  pen  movements 
(plot  size). 


Scales  data  linearly  to  fit  pen 
plotter. 


Creates  linearly  scaled  axis  with 
Hollerith  label. 


Creates  logarithmically  scaled 
axis  with  Hollerith  label. 


Draws  a svinbol  or  string  of 
SNTnbol  s . 


Draws  a straight  line  connecting  a 
set  of  linearly  scaled  points. 


DASHPT  Draws  a dashed  line  connecting  a 

set  of  linearly  scaled  points. 


Draws  a smooth  fit  connecting  a 
set  of  linearly  scaled  points. 


Draws  a straight  line  connecting 
a set  of  logarithmically  scaled 
points . 


Scales  data  logarithmically  to  fit 
pen  plotter. 


SCALOG 


Table  II 

NAME 

INIT 

TITLE 

FACTR 

NUGRID 

PLOTR 

AXES 

MAIN 

READR 


ji> . jiiMiijJ-ii'*  ii"  m.  I jmgu-.-  VM\ 


INTERPOL  Interface  Routines 


FUNCTION 


Initializes  pen  plotter 
and  sets  up  graph  origin. 


Writes  granh  title  block. 


Reduces  or  enlarges  size  of 
of  plot. 


Draws  a grid. 


Plots  all  data  and/or 
generates  histograms. 


Draws  and  labels  axes. 


Main  program  control. 


Inputs  all  data,  handles 
all  interactive  I/O. 


Checks  for  zeros  in 
logarithmically  scaled  data. 
(See  Section  III. A) 


FIXUP 


PLOTR--This  subroutine  is  the  backbone  of  INTERPOL.  All  data  is  plotted 
from  PLOTR.  An  option  available  with  INTERPOL  which  allows  for 


• histogram  generation  based  on  single  sets  of  (x,y)  data  (i.e.,  a 

new  set  of  data  is  calculated  to  enable  a histogram  st/le  plot  to 
1 be  drawn)  is  handled  via  PLOTR.  Also,  since  the  Cal  Comp  routines 

j DASHLN  and  FLINE  will  not  work  with  logarithmically  scaled  data, 

|i  PLOTR  re-scales  these  values  (if  in  logarithmic  mode)  so  as  to  be 

[{  compatible  with  the  dash  line  or  smooth  fit  options. 

All  routines  are  marked  with  comment  cards  to  aid  in  following  the 
j programming  or  instituting  changes.  A complete  FORTRAN  listing  can 

i|  be  found  in  Appendix  B. 

i' 

III.  I/O  HANDLING 

L A.  Interactive  I/O 

ij 

Subroutine  READR  handles  all  program  requests  and  user  replies. 

Each  question  is  self-explanatory  and  supplies  the  user  with  a choice  of 
r responses.  INTERPOL  precedes  each  question  with  a double  arrow(">>") 

for  clarity.  Where  a "YES"  or  "NO"  response  is  applicable,  a "Y"  or  "N" 

''  will  suffice.  Every  input  (except  for  Hollerith  data)  is  checked  for 

1 an  error,  and,  if  found,  notifies  the  user  accordingly  and  repeats  the 

question.  Warnings  or  general  information  not  requiring  user  input  are 
prececed  by  a double  asterisk  ("**"). 

:*  All  data  to  be  plotted  must  be  presented  in  (x,y)  pairs  using  free 

form  (open)  format.  While  entering  any  data  which  is  to  be  scaled 
‘ logarithmically,  INTERPOL  checks  each  (x,y)  point  for  zeros.  If  a zero 

is  encountered,  the  zero  is  reset  to  the  lowest  value  of  the  data  set 
>'  (i.e.,  the  lowest  x or  y value  input).  A warning  is  then  printed  telling 

I,  the  user  this  action  was  taken. 

I 

Several  questions  have  default  answers  and  the  user  is  supplied 
I this  information  as  necessary.  Entering  a zero,  or  transmitting  a blank 

line,  will  set  the  default  value. 

f Upon  problem  completion,  INTERPOL  will  write  the  number  of  graphs 

generated  and  compute  the  maximum  height  (y  axis  length)  encountered. 

I A message  is  printed  notifying  the  user  of  any  paper  size  restrictions. 

! A final  note  reminds  the  user  of  the  logical  unit  number  upon  which 

the  plot  records  were  written. 

' B.  Plot-Record  Tape 

t 

L INTERPOL  is  written  for  off-line  Cal  Comp  magnetic  tape  units, 

r Therefore,  the  specific  device  or  procedure  for  generating  these  tapes 

- may  vary  with  installation.  The  records  written  on  the  selected  logi- 

j cal  unit,  however,  are  all  that  are  necessary  to  drive  the  pen  plotter. 

Whether  written  directly  on  tape  or  copied  to  tape  from  mass  storage  is 
strictly  a matter  of  procedure  as  defined  by  the  user's  system.  (The 
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current  UNIVAC  1108  version  of  INTERPOL  will  dynamically  assign  a tempo-  ' 

rary  mass-storage  file  "29",  by  default,  if  so  requested.  This  may  then  I 

be  copied  to  magnetic  tape  as  mentioned  above,  and/or  copied  to  a perma- 
nent catalogued  file  upon  problem  completion).  j 

IV.  3-D  PLOTTING  PACK.AGE 

i 

A.  Tracor  Routines  i 

A thre^-dimensional  plotting  routine  developed  by  Tracor  Computing 
Corporation  was  modified  and  incorporated  for  use  with  1NTF;RP0L.  Brief-  | 

ly,  the  3-D  package  reads  sets  of  (x,y)  data  and  successively  steps  each  | 

set  in  a Z,  or  depth  direction.  Any  portion  of  data  which  falls  behind  j 

a previous  set  is  either  not  plotted,  or  plotted  with  a dashed  line  j 

(thus  giving  a "hidden"  effect).  Two  subroutines  govern  the  cal-  j 

culation  of  the  maximum  and/or  minimum  function  (i.e.,  "seen"  vs.  "unseen"  i 

data),  HIDE  and  LOOKUP.  HIDE  is  the  main  routine  and  contains  the  Cal 

Comp  software  calls.  Figure  2 is  a representation  of  the  3-D  option.  i 

B.  INTERPOL  Interface  i 

Due  to  the  special  nature  of  the  Tracor  routines  the  3-D  ca])-  > 

abilities  are  a bit  more  restrictive  than  those  previously  described, 
but  follow  the  same  format.  Only  a single  line  title  is  permitted 
(placed  above  the  plot).  Linear,  semi -log,  and  log- log  scaling  is 
available.  Only  the  x and  y axes  have  labels  and  tick  marks,  the  dejith 
axis  does  not.  There  is  no  grid  capability.  Dashed  lines  are  reserved 
for  "hidden"  data  only.  Histogram  generation  is  not  jirovided.  To 
maintain  versatility,  the  user  is  capable  of  switching  from  2-D  to  5-D 
and  back  as  successive  graphs  are  generated. 


4,  Hugli  Williamson,  "Hidden  Line  Plotting  Routine  (Algoin' thn  4P.0)", 
Cormuniaation  of  the  ACM,  Vol  IS,  No.  3 Fehnuarg  1973. 


.iJ 
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Figure  2.  Example  of  3-D  Option 
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j V . SUM^RY 

f 

! An  interactive  computer  program  for  plotting  on  off-line  Cal  Comp 

systems  has  been  developed.  The  program  acts  as  an  interface  between 
user  specified  requests  and  standard  Cal  Comp  software.  Execution  is 
j designed  primarily  for  generating  graphs  of  (x,y)  data  in  typical  two 

; dimensional  or  three  dimensional  form.  Input  is  handled  via  a question 

and  answer  conversation  with  the  user.  As  such,  an  efficient,  simple 
i and  versatile  method  of  generating  computer  graphs,  negating  the  need  for 

[ writing  special  programs  on  a case  by  case  basis,is  provided.  Little 

or  no  computer  plotting  experience  is  required  to  successfully  execute 
this  program.  Written  in  standard  FORTRAN,  this  program  is  compatible 
with  most  computing  machines. 


Tm  rzi^^urntr" 


APPENDIX  A 

Contained  in  this  appendix  are  five  plots  (Figures  A1  - AS) 
depicting  a cross-section  of  options  available  with  INTERPOL.  Also 
included  (Table  AI)  is  a list  of  integer  equivalents  of  Cal  Comp 
symbols • 


FIGURE  fl- 


X AXIS  LABEL 

Figure  A1 . Smooth  Fit,  Log-Log  Scaling, 
Vertical  Grid  Lines 


Figure  A3.  Histogram,  Linear  X-Log  Y Scaling 


FIGURE  fl- 


Figure  A4.  Points  Connected  by  St 
Linear  X and  Y Scaling 


Figure  A5.  Dashed  Line,  Log-Log  Scaling 
Vertical  and  Horizontal  Grid 
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Table  AI.  Integer  Equivalents  of  Cal  Comp  Symbols 


C I*W*T*E*R*P*0*L 

C 

C A GENERAL  PURPOSE  2-DIf-BiNSIONAL/3-DINENSIONAL  PLOTTING  ROUTINE 

C 

C 

C 

C 

C THIS  IS  TIIE  UNIVAC  1108  VERSION. 

C 

C MAIN  ROUTINE  FOR  'INTERPOL': 

C 

C IHE  ARRAYS  FOR  THE  (X.Y)  DATA  ARE  ARBITRARILY  SET  TO  900 

C IN  COMMON  DATA.  THE  NUMBER  OF  SETS  OF  DATA  PER  GRAPH 

C IS  ARBITRARILY  SET  TO  20  VIA  COMMON  DRAIV. 

C THE  DIMENSION  STATENENT  (ICHK)  , BELOW  , WAS  ARBITRARILY  SET 

C TO  100  --  I.E.,  AN  ARBITRARY  NUMBER  OF  GRAPHS  PER  RUN. 

C IT  IS  SIMPLY  USED  TO  QIECK  VERTICAL  AXIS  LENGTHS  FOR  PAPER 

C SIZE  RESTRICTIONS. 

C 

C...  THE  NEXT  CARD  IS  FOR  CDC  MACHINES. 

C 

C PROGRAM  MAIN(INPUT, OUTPUT, TAPE5=INPUT,TAPE6=OUTPUT,TAPE29) 

C 

CONWON  /INPUT/  IXAXIS(5),IYAXIS(5),ITITLI(5),ITITL2(5),ITITL3(5), 
1ITITL4(5),ITITL5(5),ITYPE,XAXIS,YAXIS,FACT,ITITLE,Y2,X1,IGRID, 
2YINC,XINC,N0U,NIN,ILINES,LDEV,NSETS,IHN,IHY,IHB,IH0 
C 

CONMON  /DATA/  XARRAY (900) ,YARRAY (900) ,HIST0X(900) ,HIST0Y (900) , 
IXPTS  (900)  ,YPTS  (900)  ,ANUX(900)  ,ANUY  (900)  ,A;>iUHX(900)  .ANUHY  (900) 

C 

COMMON  /DRAW/  NPTS (20) , I PLOT(20) , IHIST (20) , LINTYP (20) , INTEQ (20) , 
1LOGTYP(20) .SKIP 
C 

DATA  NOU,NIN/o,5/ 

DATA  IHN,IHY,IHB,  1110/ IHN.IHY,  111  ,1110/ 

C 

DIMENSION  laiK(lOO) 

WRITE  (NOU,5000) 

SKIP=0. 

NPASS=1 

100  CALL  READR  (NPASS.MASK) 

CALL  AXES  (MASK) 

IF  (MASK.EQ.IHY)  GO  TO  105 
IF  (IGRID.EQ.IIIY)  CALL  NUGRID 
IF  (ITITLE.EQ.IHY)  CALL  TITLE 
CALL  PLOTR 
GO  TO  no 

lOb  CALL  HIDE  (YSTEP) 

YAXIS=YAXIS+YSTEP 
no  iaiK(NPASS)=YAXIS*FACT 


28 


n n n n n n n n n n n n n 


11811  WRITTEN  ON  UNIT  ,12.211  .) 

5040  FORMAT  (/30H***  PLOT(S)  WILL  FIT  ON  NARROW  PAPER.) 

5050  FORMAT  (/4311***  PLOT(S)  WILL  NOT  FIT  ON  NARROW 

PAPER,, /,4X, 

IlbHUSE  WIDE  SIZE.) 

5000  FORMAT  (/2011***  END  OF  RUN /) 

END 

SUBROUTINE  AXES  (MASK) 

THIS  SUBROUTINE  SCALES  ALL  DATA  AND  DRAWS  THE  APPROPRIATE  AXES 
WITH  LABELS  AS  WELL  AS  THE  GRAPH'S  BORDER. 


CONWON  /INPUT/  IXAXIS(5),IYAXIS(5),ITITLI(5),ITITL2(5),ITITL3(5), 
1ITITL4(5),ITITL5(5) ,ITYPE,XAXIS,YAXIS,FACT,ITITLE,Y2,XI,IGRID, 
2YINC,XINC,N0U,NIN,ILINES,LDEV,NSETS,IHN,IHY,IHB,IH0 

COMN»N  /DATA/  XARRAY (900) .YARRAY (900) ,HISTOX(900) .HISTOY (900) , 
IXPTS(900)  ,YPTS  (900)  ,ANUX(900) , ANUY  (900) , ANU!IX(900)  .ANUllY  (900) 

CONWON  /DRAW/  NITS  (20) , I PL0T(20) , I HIST  (20) , LIN  IT'P  (20) , INTEQ(2U) . 
1LOGTYP(20) .SKIP 

ITOT=0 

DO  100  I=I.NSETS 
ITOT=ITar+NPTS(I) 

100  CONTINUE 

GO  TO  (105. no. 115. 120) . ITYPE 

...  ITYPE  = 1 LINEAR  X.LINE/\R  Y SCAUil)  A.XES. 

105  CALL  SCALE  (XARRAY . XAXIS . ITOT, 1) 

CALL  SCALE  (YARRAY.YAXIS. ITOT. I) 

IF  (MASK.EQ.IHY)  GO  TO  125 

CALL  AXIS  (0. O.O.O.IXAXIS. -25. XAXIS. 0. 0. XARRAY(IT0T+1) . 
IXARRAY(IT0T+2)) 

CALL  AXIS  (0.0.0.0.IYAXIS.25.YAXIS.90.0.YARRAY(IT0T+1) . 
IYARRAY(IT0T+2)) 

CALL  PLOT  (0.0.YAXIS.3) 

CALL  PLOT  (XAXIS.YAXIS.2) 

CALL  PLOT  (XAXIS. 0.0. 2) 

GO  TO  125 

...  ITYPE=2  LINEAR  X.LOG  Y SCAUID  AXES. 

no  CALL  SCALE  (XARRAY  .XAXIS.  ITOT,  I) 

CALL  SCALG  (YARRAY.YAXIS, ITOT, I) 

IF  (MASK.EQ.IHY)  GO  TO  125 

CALL  AXIS  (0. 0.O.O.IXAXIS, -25, XAXIS, 0. 0, X.4RRAY(IT0T+I)  , 

29 


n n n n non  non 


■JiUlJ  -JlUJPi 


i 

i 


1 


i 


il 


I 


T 


iXARRAY(IT0T*2}) 

CALL  LGAXS  (O.O.O.O, IYAXIS,25,YAXIS,90.0,YARRAY (ITOT+1) , 
lYARRAY(IT0T+2)) 

CALL  PLOT  (0.0,YAXIS,3) 

CALL  PLOT  (XAXIS,YAXIS,2) 

CALL  PLOT  (XAX1S,0.0,2) 

GO  TO  125 

. . ITYPE=3  LOG  X, LINEAR  Y SCALED  AXES. 

15  CALL  SCALG  (XARRAY.XAXIS, ITOT, 1) 

CALL  SCALE  (YARRAY,YAXIS,1T0T, 1) 

IF  (MASK.EQ.IHY)  GO  TO  125 

CALL  LGAXS  (0.0 ,0 . 0 , I XAXIS, -25 , XAXIS , 0 . 0, XARRAY (ITOT+ 1) , 
1XARRAY(1T0T*2)) 

CALL  AXIS  (0.0,0.0,IYAXIS,25,YAXIS,90.(),YARRAY(ITOT+1)  , 
lYARRAY(ITOT+2)) 

CALL  PLOT  (0.0,YAXIS,3) 

CALL  PLOT  (XAXIS,YAXIS,2) 

CALL  PLOT  (XAXIS,0.0,2) 

GO  TO  125 


...  ITYPE  = 4 LOG  X.LOG  Y SCALED  AXI-S. 

120  CALL  SCALG  (XARRAY, XAXIS, ITOT, 1) 

CALL  SCALG  (YARRAY.YAXIS, ITOT, 1) 

IF  (MASK.EQ.IHY)  GO  TO  125 

CALL  LGAXS  (0 .0, 0 . 0, I XAXIS, -25 , XAXIS , 0. 0, XARRAY (ITOT+ 1) , 
lXARRAY(ITOT+2)) 

CALL  LGAXS  (0 . 0 ,0 . 0 , lYAXIS , 25 ,YAXIS ,90. 0, YARRAY (ITOT^ 1) , 
lYARRAY(ITOT+2)) 

CALL  PLOT  (0.0,YAXIS,3) 

CALL  PLOT  CXAXIS,YAXIS,2) 

CALL  PLOT  (XAXIS, 0.0, 2) 

125  RETURN 
END 

SUBROUTINE  FACTR 

THIS  SUBROUTINE  CHANGES  THE  LENGTH  OF  ALL  PEN  MOVEMENTS  BY  A 
USER  SUPPLIED  SCALING  FACTOR  (SUBR.  RE,\DR)  . 1.0  IS  FULL  SCALE. 

COMMON  /INPUT/  I XAXIS (5) , I YAXIS (5) , ITITLl (5) , 1TITL2 (5) , 
1ITITL3(5),ITITL4(S3 ,ITITL5(5) ,ITYPE,XAXIS,YAXIS,FACT,ITITLE,Y2,X1, 
2IGRIU,YINC,XINC,NOU,NIN,  ILINES,  LDI-V,NSETS,  IIW,IHY,  IHB,  1110 
C 

COMMON  /DATA/  XARRAY  (900)  ,YARRAY  (900)  ,HIST0X(900)  ,1IIST0Y  (900)  , 
1XPTS(900) ,YPTS (900) ,ANUX(900) ,ANUY (900) ,ANU11X(900) ,ANUHY (900) 

C 

COMMON  /DRAW/  NPTS  (20) , IPLOT(20) , IHIST(20) , LINTYP(20)  , II'ITEQ(20)  , 
ILOGTYP(20),SKIP 


n n n n n n r;  n n.  n n 
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CALL  FACTOR  (FACT) 

RETURN 

END 

SUBROUTINE  FIXUP 

THIS  ROUTINE  SETS  ANY  DATA  WHICH  IS  TO  Bi;  LOCAR ITHMICALLY  SCALED 
AND  WAS  FOUND  TO  BE  EQUAL  TO  ZERO  TO  THE  MINIMUM  OF  THE  RESPECTIVE 
X OR  Y ARRAY.  THIS  PREVENTS  AN  ABORTED  RUN.  THE  USER  IS  TOU). 

COMMON  /INPUT/  IXAXIS{5) , IYAXIS(5) , ITITLI(5) , ITITLJ(5) ,rriTL3(5) , 
IITITL4(5),ITITL5(5) ,ITYPE,XAXIS,YAXIS,FACT,ITITLE,Y2,X1,IGRID, 
2YINC,XINC,NOU,NIN,ILINES,LDEV,NSETS,IHN,IHY,IHB,IliO 

COMMON  /DATA/  XARRAY(900)  .YARRAY  (900)  ,11IST0X(900)  .HISTOY  (900) , 
IXPTS  (900)  ,YPTS  (900),  ANUX(900),  ANDY  (900)  ,ANUiIX(900)  ,ANUHY(900) 

COMMON  /DRAW/  NPTS (20) , IPLOT(20) , IHIST (20) , LI NTYP (20) , lN-n;Q(20) , 
1LOGTYP(20),SKIP 

DO  1 I=I,NSETS 
ITOT=rrOT+NPTS(I) 

1 CONTINUE 
XMIN=1.0E+29 
YMIN=1.0E+29 

DO  2 I=I,ITOT 
XMIN=AMIN1(XMIN,XARRAY(1)) 

YMIN=AM1N 1 (YMIN .YARRAY (1) ) 

2 CONTINUE 

IF  (XMIN.EQ.O.)  XMIN=XMIN+.0001 

IF  (YMIN.EQ.O.)  YMIN=YMIN+.0001 

L=rrYPE-I 
GO  TO  (3,S,7) , L 

j DO  4 I=I,ITOT 

IF  (YARRAY (I ).EQ.O.)  YARRAY ( I ) =YMIN 
J CONTINUl; 

GO  TO  9 

a DO  0 i=i,rroT 

IF  (XARRAY(I) .EQ.O.)  XARRAY(1 )=XMIN 
t)  CONTINUl- 

GO  TO  9 

7 DO  8 I=l,ITOT 

IF  (URRAY(I).EQ.O.)  XARRAY  ( I ) =XMIN 
IF  (YARRAY(I).EQ.O.)  YARRAY (I )=YMIN 
< CONTINUE 

0 CONTINUE 

WRITE  (NOU, 10) 

RliITJRN 


*****FORMAT  stati;mi;nts***** 
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10  FORMAT  (/Sill***  INPUT  WARNING:  ZHRO'S  LNCOUNniRRl)  ON  LOGARITUMI . 
119I1CALLY  SCALRU  DATA- - , / , 3ol  1 VALURS  WRRR  RR-SRT  TO  MINIMUM  OF, 

23311  RRSPECTIVR  ARRAY  (XMIN  OR  YMIN).) 

RNU 

SUBROUTINR  IIIDR  (YSTRP) 

nils  ROUTINR  PLOTS  niR  3-l)IMRNSIONAL  GRAPH. 

COMMON  /INPUT/  IXAXIS(5) .lYAXISCS) ,ITITL1(5) , ITITL2(5) , ITITL3(rO  , 
1ITITL4(5),ITITL5C5),ITYPR,XAXIS,YAXIS,FACT,ITITLR,Y2,X1,IGRI1), 
2YINC,XINC,NOU,NIN,ILINES,Ll)RV,NSRTS,IIlN,im',IIIB,  IMO 

COMMON  /DATA/  XARRAY(900)  .YARRAY  (900)  ,I1ISTOX(900)  .IIISTOY  (900) , 
IXPTS  (900)  ,YPTS  (900) , ANUX(900)  ,ANUY  (900)  ,ANII1X(900)  .ANUIIY  (900) 

COMMON  /ORAlV'/  NPTS  (20) , IPLOT(20) , I1IIST(20)  ,LINTYP(20)  , INTI;Q(-!0)  , 
1L0GTYP(20),SKIP 

OIMRNSION  X(900)  ,Y  (900)  , XG(900) , G(900) , XII  (900)  ,11(900) 

RQUIVAIiiNCR  (Kl.IWIlICIl),  (K2,SLOrR),  (FNSMl.Zl),  (IGGP1,K1),  (Kl, 
1N2) 

DATA  RPS1,MAXUIM,XSTART,YSTART/1 .R-9,900,4. , 3. / 
F(XX,XI,YI,XIP1,YIP1)=YI+(XX-XI)*(YIP1-YI)/(X1P1-XI) 

WRITR  (NOU.SOOO) 

RRAD  (NIN,5  0 30)  I DASH 

YSTRP=YSTART 

SKIP=XAXIS+4.0 

JN=1 

J0=0 

NG=0 

rroT=o 

DO  100  I=1,NSRTS 
rrOT=ITOT+NPTS(l) 

100  CONTINUR 
C 

DO  355  JM=1,NSRTS 
N1=NPTS (JM) 

JO=JO+Nl 

IK=1 

IX)  130  IJ»JN,JO 
rrOT=ITOT 

GO  TO  (120,  105,110,  115),  m'PF. 

105  Y (IK)  = (L0G10(YARRAY  (1.1)  )-L0G10(YARRAY  (rrOT+  1) )) /YARRAY  ( rrOT+ 

1 ) 

(IK)-XARRAY(I.I) 

GO  TO  125 

1 1 0 X ( I K)  = ( LOG  10  (XARRAY  ( 1.1) ) - U)G  10  ( XARR AY  ( 1 TOl>  1 ) ) ) /X ARRAY  ( 1 T01> 

1 2) 

32 


I 


i 

i 


1 

j 


Y(IK)»YARRAYC1J) 

GO  TO  125 

115  X{IK)  - (LOGIO  (XARRAY  ( I J ) ) - LOG  10  (XARRAY  ( ITOT-^  1) ) ) /XARRAY  (1X01  + 

1 2) 

Y ( I K ) - 1 LOG  1 0 C Y ARR AY  ( IJ ) ) - me  1 0 (Y ARRAY  ( 1 TOT+ 1) ) ) /Y ARR AY  ( I TOT+ 
1 2) 

GO  TO  125 

120  X (IK) -XARRAY (IJ) 

Y(IK)-YARRAY(IJ) 

125  IK-IK+1 

130  CONTINUE 

JN-JN+Nl 

DO  135  1-2, N1 

IF  (X(I-l) .LT.X(I))  GO  TO  135 
WRITli  (NOU,5010)  X{I-1),X(I) 

RLTURN 

135  CONTINUi: 

IF  (JM.GT.l)  GO  TO  205 

NFNS-NSFTS 

XMIN-XARRAY(1T0T+1) 

YMIN-YARRAY(1T0T+ 1) 

OF  LTAX- XARRAY  ( 1X01'+  2 ) 

01:LTAY-Y ARRAY  (1X01'+  2) 

IF  CNU4.  U;.MAX1)1M)  GO  TO  140 
GO  TO  30 0 
140  SIGN-1. 

IF  (NG.LT.-l)  SHa^  = -l. 

IF  (NG.FQ.-1.OR.NG.F0.-3)  GO  TO  145 
CALL  PLOT  (0.,YST/\RT+YAX1S,3) 

CALL  OASllP  (XAXIS,YSTART-t-YAXIS,.03) 

CALL  DASllP  (XAXIS,YSTART,.03) 

CALL  OASllP  (0.  ,YSTART,.03) 

CALL  PLOT  (XAXIS,YSTART,3) 

CALL  OASllP  (XAXIS-fXSTART.0.,.03) 

145  CjM.L  symbol  (2. 0,YAX1S+YSTART+.  10,.  14,111X1.1,0. ,25) 

GO  TO  (150,155,100,  105),  ITYPi: 

150  i:ALL  axis  (XST.\RT,0.  ,IXAX1S,-25,XAX1S,0.  ,XMlN,Oi;i,TAX) 

CALL  PLOT  (.XSTART,0.  ,3) 

CALL  PLOT  (0.,YSTART,2) 

CALL  AXIS  (0.,YSTART,IYAX1S,25,YAXIS,00.,YMIN,0I:LTAY) 

TO  170 

155  CALL  AXIS  (XSTART,0 . , I XAX1S,-25,XAXIS  ,0 . , XM1N,01.LTAX) 

CALL  PLOT  (XSTART,0.,3) 

CALL  PLOT  (0.,YSTART,2) 

CALL  LG>\XS  (0.  ,YSTART,1YAX1S,25,YAXIS,'.)0.  ,YM1  N ,01;LTAY  ) 

GO  TO  170 

loO  CALL  LGAXS  (.XSTART, 0 . . 1 XAXl S , -25,  XAXl S, 0 . , XMIN , OFLTAX) 

CALL  PLOT  (XSTART,0.,3) 

CALI.  PLOT  (0.  ,YSTART,2) 

CALL  AXIS  ( 0 . , Y START , 1 YA X I S , 25  , Y AX  1 S , 0 0 . , Y Ml  N , Oil LTAY ) 

U)  TO  170 
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IbS  CALL  LGAXS  (XSTART.O. , IXAXIS,-25,XAXIS,0 . .XMIN.DELTAX) 
CALL  PLOT  9xstart,0.,3) 

CALL  PLOT  (0.,YSTART,2) 

CALL  LGAXS  (0. , YSTART, I YAXIS,25 , YAXIS,90. , YMIN.DELTAY) 
170  INDEXT-3 

GO  TO  (190,175,180,185),  ITYPE 
175  YMIN-0. 

OELTAY-1. 

GO  TO  190 
180  XMIN=«0. 

UELTAX-1. 

GO  TO  190 
185  XMIN-0. 

l)ELTAX»l. 

YMIN-0. 

UELTAY-1. 

190  CONTINUE 

IF  (NFNS.LE.O)  GO  TO  195 
FNSMl-NFNS-1 

0XIN»XSTART*DELTAX/FNSM1 
l)YIN=YSTART*DELTAY/FNSMl 
195  00  200  J=1,N1 

XG (INDEX!) =X(J) 

G(1NUEXT)=SIGN*Y(J) 

INUEXT=INDEXT+1 
200  CONTINUE 

EPS=EPSl*(ABS(XMlN)-*-ABS(DELTAX) ) 

NG-NU4 

XG(1)  = -FNSM1*DXIN+XMIN-ABS (XMIN)-ABS(Xn(3) )-l . 

XG(2)=XG(3)-EPS 

XG  (Nl-t-3)=XG(Nl  + 2)+EPS 

ZZ-YMIN 

IF  (SIGN.LT.O.)  ZZ--YMIN-50.*DELTAY 

G(1)-ZZ 

G(2)=ZZ 

a(Nl+3)-ZZ 

G(NG)»ZZ 

CALL  PLOT  (XSTART,0,,-3) 

X(Nl-t'l)-XMIN 

X(Nl-*-2)-UELTAX 

Y(N1+1)«YMIN 

Y(NI+2)-UELTAY 

CALL  LINE  (X,Y,N1, 1,0,0) 

UXKK-0. 

DYKK-0. 

RELINC-DliLTAX/DELTAY 

XG(NG)-SIGN 

GO  TO  355 

205  SIGN-XG(NG) 
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XGCNG)=X(N1) 

IF  (NFNS)  225,215,210 
210  UXKK=»UXKK+DXIN 
DYKK=l)YKK+DYIN 
215  DO  220  J=1,N1 

Y(J)=»SIGN*(Y(J)+DYKK) 

X(J)=X(J)-DXKK 
220  CONTINUE 

225  CALL  LOOKUP  CXU),XG(1)  ,JJ) 

IF  (JJ.GL.MAXDIM)  GO  TO  3oO 
DO  230  J=1,JJ 
XH(J)=«XG(J) 
ll(J)=G(J) 

230  CONTINUL 

IG=JJ+1 
X11CIG)=X(1) 

H(IG)=F(X(1),XG{JJ)  ,G(JJ),XG(IG),G(IG)) 

INDl-XG=JJ 

INDFXT=1 

Z1=X(1) 

Fl=ll(IG)-Y(l) 

IT=2 

JJ=1G 

IF  (}1(IG)  .GF.Y(l)  ) GO  TO  235 
IF  (JJ.Gli.NlAXDIM)  GO  TO  3i)0 
JJ=IG+1 
I1CJJ)=Y(U 
Xll(JJ)=ZUliPS 
235  1AST=0 

X1  = Z1 

240  IF  (XG(IG)  .LT.X(rr))  GO  TO  245 

IWI1ICU=0 

X2=x(rr) 

F2=F(X2,XG(IG-1),G(IG-1),XG(IG),G(1G))-Y(1T) 
IT=lT-t-l 
GO  TO  250 
245  X2=XG(IG) 

IW11ICH=1 

F2=G(IG)-F(X2,X(rr-l),YHT-l),X(lT),Y(rn) 

IG=IG+1 

250  IF  (FI*F2.GT.O.)  GO  TO  200 

di;nom=x2-xi 

IF  (DliNOM.FQ.O.)  I)I;N0M=  . 00001 
SL0PF=(F2-F1)/D1'.N0M 
IGG=IG-1-IWIIIC11 
ITT=IT-2+IWIliai 

IF  (AUS(SLOPF.*IU'.LINC)  .GT.  GO  TO  255 

Z2=X2 

GO  TO  270 


3S 
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260 

265 

270 


275 

280 

285 


290 


Z2=X1-F 1/SLOPE 
GO  TO  270 
X1=X2 
F1=F2 

IF  (IT.LE.Nl)  GO  TO  240 

IAST=1 

Z2=X  (Nl) 

CALL  LOOKUP  CZ2,XG(INDEXG) , IGG) 

IGG=INDEXG+IGG-1 

ITT=N1-1 

ZZ=.99*Z1+.01*Z2 

CALL  LOOKUP  (ZZ,X(INDEXT) ,K1) 

CALL  LOOKUP  (ZZ,XG(INDEXG) ,K2) 

Kl»Kl-*-INDEXT-l 

K2=K2+INDEXG-1 

IF  (F(ZZ,X(K1) ,Y(K1) ,X(K1+1) ,Y (Kl+1) ) .GT.F(ZZ,XG(K2) ,G(K2) , 

XG(K2+1) ,G(K2+1)))  GO  TO  300 

IF  (JJ+IGG-INDEXG.GE.MAXDIM)  GO  TO  360 

NGR=ITT-INDEXT+2 

NN2=JJ 

NJJ=JJ 

ANUX(NJJ)=XH(NJJ) 

ANUY  (NJJ)  =H  (NJJ) 

IF  (NGR.EQ.2)  GO  TO  280 
NJ1=INDEXT+1 

DO  275  I=NJI,ITT 
NJJ=NJJ+ I 
ANUX(NJJ)=X(I) 

NUY{NJJ)=Y(I) 

CONTINUE 

NJJ=NJJ+I 

ANUX(NJJ)=Z2 

ANUY (NJJ) =F(Z2,X (ITT) ,Y (ITT) ,X(ITT+1) , Y (ITT+1) ) 

NM=NN2+NGR-I 

NL=0 

DO  285  I=NN2,NNM 
NL=NL+I 

ANUHX(NL)=ANUX(I) 

ANUUY(NL)=ANUY(I) 

CONTINUE 

ANUHX(NL+1)=XMIN 

ANUHY(NL+1)=YMIN 

ANUHY(NL+2)=DELTAY 

ANUHX(NL+2)=DELTAX 

IF  (IDASH.EQ.IHY)  CALL  DASHL  (ANUHX,ANUUY,NL, I) 

IF  (INDEXG.EQ.IGG)  GO  TO  295 
J1=INDEXG+1 

DO  290  I»J1,IGG 

JJ=JJ+1 

XH(JJ)=XG(I) 

U(JJ)=G(I) 

CONTINUE 


I 
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295  JJ=JJ+1 

XH(JJ)=22 

II(JJ)=F(Z2,XG(IGG),G(IGG),XG(IGG+1),G(1GG>1) ) 

INDFXG-IGG 

INDEXT»ITT 

GO  TO  320 

300  NGRAPH=ITT-INUEXT-t-2 

IF  (JJ+NGRAPH-l.GT.MAXDlM)  GO  TO  300 
N2=JJ 

IF  (NGRAPH.EQ.2)  GO  TO  310 
J1=*INDEXT+1 

DO  305  I=J1,ITT 

JJ=JJ+1 

XH{JJ)=X(I) 

H(JJ)=Y(I) 

305  CONTINUE 

310  JJ=JJ*1 

X11(JJ)=Z2 

UCJJ)=F(Z2,X(ITT),Y  (ITT)  .XCITT^D.YCITT-^l)) 

NM=N2+NGRAPH-1 

L=0 

DO  315  I=N2,NM 
L=L+1 

XPTS(L)=X)J(I) 

YPTS(L)=U(I) 

315  CONTINUE 

XPTS  (L>1)=XMIN 

XPTS(L+2)=UELTAX 

YPTS(L+1)=SIGN*YMIN 

YPTS(L+2)=SIGN*DELTAY 

CALL  LINE  (XPTS, YPTS.L, 1,0,0) 

INDEXT=>ITT 

INDEXG=IGG 

320  IF  (LAST.EQ. 1)  GO  TO  325 

X1=X2 
F1=F2 
Z1  = Z2 

IF  (IT.LE.Nl)  GO  TO  240 
GO  TO  2t>5 

325  IF  (XG(NG) .LE.XG(NG-l))  NG=NG-1 

IF  (XG(NG) .LE.X(Nl))  GO  TO  335 
IF  (JJ+3+NG-IGG.GT.MAXDIM)  GO  TO  3o() 
XH(JJi-l)»Xll(JJ)+EPS 
JJ= JJ+ 1 

H ( J J)  =F  (X  (N 1) , XG  ( I GG) , G ( I GG)  , XG  ( 1 GG*  1 ) , G ( I GG+ 1 ) ) 
IGGPI=IGG+1 

DO  330  J=1GGP1,NG 
JJ-JJ+1 
XH(JJ)*XG(J) 
ll(JJ)»G(J) 
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330  CONTINUE 

335  NG-JJ+2 

IF  (NG.GE.MAXDIM)  GO  TO  3b0 
DO  340  I»I,JJ 
G(I)-H(I) 

XG(I)=XU(I) 

340  CONTINUE 

XG(JJ+1)*XG(JJ)+EPS 

G(JJ+I)»YMIN+DYKK 

IF  (SIGN.LT.O.)  G(JJ-^I)=-YMIN-50.*DELTAY+DYKK 
G(NG)-G(JJ+1) 

IF  (NFNS.LT.O)  GO  TO  350 
DO  345  1=1, N1 
X(I)=X(I)+DXKK 
Y(I)=SIGN*Y(I)-DYKK 
345  CONTINUE 

350  XG(NG)=SIGN 

355  CONTINUE 

CALL  PLOT  (XAXIS,0.0,3) 

CALL  PLOT  (SKIP, 0.0, -3) 

RETURN 

3t)0  WRITE  (6,5020)  MAXDIM 
ETURN 


**** ‘FORMAT  STATEMENTS***** 


5000  FORMAT  (/37H  DO  YOU  WANT  ANY 


HIDDEN 


-D  data,/4X, 

140HPL0TTED  WITH  A DASHED  LINE  --  YES  OR  NO?) 

5010  FORMAT  (/33H***  INPUT  ERROR  (HIDE):  X(I-l)  = , IPEIO .4, 5X, 6HX (I ) =, 
IIX, 1PE10.4,/,35H  X(I-l)  MUST  BE  LESS  THAN  X(I)) 

5020  FORMAT  (/48H***  ERROR  IN  DIMENSIONED  ARRAYS  (HIDE):  MAXDIM  =,I4,/, 
135H  INCREASE  MAXDIM  TO  RUN  PROBLEM.) 

5030  FORMAT  (Al) 


END 

SUBROUTINE  INIT 


THIS  SUBROm'INE  INITIALIZES  THE  PLOTTER.  THE  (X,Y)  COORDINATES 
IN  THE  DATA  STATENENT  SET  THE  POSITION  OF  THE  GRAPH'S  ORIGIN 
ON  THE  PAPER  (I.E.,  LOWER  LEFT  HAND  CORNER). 

COMMON  /INPUT/  I XAXIS (5) , lYAXIS (5) , ITITLl (5) , ITITL2 (5) , ITITL3 (5) , 
1ITITL4(5) ,ITITL5(5) ,ITYPE,XAXIS,YAXIS,FACT,ITITLE,Y2,X1,IGRID, 
2YINC,XINC,NOU,NIN,ILINES,LDEV,NSETS,IHN,IHY,IHB,IHO 
C 

COMMON  /DATA/  XARRAY  (i)00)  ,YARRAY  (900)  ,HISTOX  (900)  ,HISTOY  (900) , 
1XPTS(900),YPTS(900),ANUX(900),ANUY(900),ANUHX(900) ,ANUHY (900) 
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common  /DRAW/  NPTS  (20) , I PLOT{20) , II1IST{20) , LINTyp(20) , INTliQ  (20)  , 
1LOGTYP(20) .SKIP 
C 

DATA  X,Y/3. 0,2.0/ 

C 

CALL  PLOTS  (IBUFF.NLOC.LUEV) 

CALL  PLOT  (0.0, -36. 0,-3) 

CALL  PLOT  (X,Y,-3) 

RETURN 
END 

SUBROUTINE  LOOKUP  (X,XTBL,J) 

THIS  SUBROUTINE  IS  CALLED  FROM  HIDE  AS  AN  EFFICIENT  METHOD 
OF  TABLE  LOOKUP  (STOLEN  FROM  TRACOR  3-D  PLOT  PACKAGE). 

DIMENSION  XTBL(I) 

J=2 

1 IF  (XTBL(J)-X)  2,3,4 

2 J=J+I 
GO  TO  I 

3 RETURN 

4 J=J-I 
RETURN 
END 

SUBROUTINE  NUGRID 
C 

C THIS  SUBROUTINE  DRAWS  GRID  LINES  AS  REQUESTED  BY  THE  USER.  IT 
C WILL  NOT  DRAW  THROUGH  THE  TITLE.  ZINC  IS  THE  INCREMENT  BY  WHICH 
C THE  PLOT  PEN  CHECKS  TO  SEE  IF  IT  IS  WITHIN  THE  TITLE  BOUNDARIES 
C SET  UP  BY  SUBROUTINE  TITLE. 

C 

COMMON  /INPUT/  IXAXIS (5) , lYAXIS (5) , ITITLl (5) , ITITL2 (5) , ITITL3(5) , 
1ITITL4(5),ITITL5(5),ITYPE,XAXIS,YAXIS,FACT,ITITLE,Y2,XI,IGRID, 
2YINC,XINC,N0U,NIN,ILINES,LI)EV,NSETS,IHN,IHY,IHB,IH0 
C 

COMMON  /DATA/  XARRAY (900) , YARRAY (900) .HISTOX (DOO) .HISTOY (900) , 
IXPTS  (900)  ,YPTS(900)  ,ANUX(900)  ,ANUY  (900)  ,ANUHX(900)  ,A.NUHY  (900) 

C 

CONWON  /DRAW/  NPTS  (20) , IPLOT(20)  , IHIST(20)  , LINTYP(20)  , INTT.Q  (20) , 
1LOGTYP(20),SKIP 
C 

DATA  ZINC/0. 1/ 

C 

IF  (ILINES.EQ.O)  GO  TO  125 
GO  TO  (100,105,110,115,120),  ILINES 
100  Yl=Y2-.3 

GO  TO  125 
105  Yl=Y2-.6 

GO  TO  125 
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110  Yl-YJ-.y 
GO  TO  125 
115  Y1»Y2-1.2 

GO  TO  125 
120  Y1-Y2-1.S 

125  CONTINUE 
X2-XW3.75 
C 

C...  DRAWS  VERTICAL  GRID  LINES. 

C 

ITOT*0 

DO  130  I-1,NSETS 
ITOT»nOT>NPTS(I) 

130  CONTINUE 

IF  (XINC.EQ.0.0)  GO  TO  loO 
IF  (ITYPE.EQ.l.OR.  ITYPi;.EQ.2)  XINC=1. 

IF  (ITYPE.EQ.3.0R.ITYPE.EQ.4)  XINC=1 . /XARRAY (ITnT+2) 
X=XINC 

ICOUNT«XAXIS/XINC 
JCOUNT=YAXIS/ZINC 
DO  155  I=«1,IC0UNT 
Y»0.0 

CALL  PLOT  (X,Y,3) 


DO  150  J=1,JC0UNT 
Y=Y>ZINC 

IF  (X.GE.X1.AND.X.LE.X2) 
GO  TO  140 

GO 

TO  13o 

135 

IF  CY.GE.Y1.AND.Y.LE.Y2) 

GO 

TO  145 

140 

CALL  PLOT  (X,Y,2) 

GO  TO  150 

145 

CALL  PLOT  (X,Y,3) 

150 

CONTINUE 

X*X>XINC 

155 

CONTINUE 

...  DRAWS  HORIZONTAL  GRID  LINES. 

IF  (YINC.EQ.0.0)  GO  TO  190 
100  IF  ClTYPE.EQ.l.OR.ITYPE.EQ.S)  YINC=1. 

IF  (ITYPE.EQ.2.0R.ITYPE.EQ.4)  YINC=1 . /YARRAY (ITOT+2) 
Y-YINC 

KCOUNT-YAXIS/YINC 
,,  LCOUNT-XAXIS/ZINC 

1)  DO  185  K-1,KC0UNT 

''  X-0.0 

ji  CALL  PLOT  (X,Y,3) 

;1  DO  180  L«1,LC0UNT 

^ X-X+ZINC 

IF(Y.GE.Y1.AND.Y.LE.Y2)  GO  TO  l&O 
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GO  TO  170 

lo5  IF  (X.GH.X1.ANI).X.U;.X2)  GO  TO  175 

170  CALL  PLOT  (X,Y,2) 

GO  TO  180 

175  CALL  PLOT  tX,Y,3) 

180  CONTINUL 

Y=Y+YINC 
185  CONTINUL 
100  RLTURN 
UNO 

SUBROUTINL  PLOTR 
C 

C THIS  SUBROUTINL  SliTS  UP  TllL  DATA  ARRAYS  WITH  APPROPRIATH  SCALING 
C FACTORS  tCOMPUTLD  IN  SUBR.  AXILS)  .GliNFRATIiS  HISTOGRAMS  IF  RLQUI'.STliD 
C AND  TURN  PLOTS  ALL  DATA  AS  INPUT  IN  SUBROUTINIi  RFADR. 

C 

C THIS  ROUTINE  IS  ALSO  NEXT  TO  INn’OSSIBLF.  TO  EXPLAIN  OR  FOLLOW 
C 

COMMON  /INPUT/  I XAXI S (5) , I YAXI S CS) , ITITLl (5) . ITITL2 (5) , ITITLS (5) , 
1ITITL4(5) , ITrrL5(5)  ,ITYPF.,XAXIS,YAXIS,FACT,rnTLE,Y2,Xl,IGRID, 

2Y I NC, XINC , NOU , N I N , I L I NFS , LDEV , NSETS , I HN , 1 1 lY , I HB , I HO 


COMMON  /DATA/  XARRAY  (000)  ,Y/\RRAY  (1)00) , 11 1 STOX  (DOO)  .HISTOY  (1)00) , 
IXPTS  (900),  YFl'S  (1)00),  ANUX(l)OO),  ANDY  (1)00)  ,ANU11X(900)  .ANU  lY  (DOO) 

COMMON  /DRAW/  NPTS (20) , irLOT( 20) . IHIST(20 ) . LINTYP (20) , INTEQ (2t) ) , 
1L0GTYP(20),SKIP 


SKIP=XAXIS+4,0 
JN  = 1 
J0=0 
rroT=o 

DO  100  1=1, NSETS 

ITOT=ITOT+NPTS(I) 

CONTINUE 


DO  255  JM=1. NSETS 
L=IPLOT(JM) 

JO=JO+NPTS (JM) 

IK=1 

DO  130  IJ=JN,.IO 
XPTS(IK)=XARRAY(IJ) 

YPTS(IK)=YARRAY(IJ) 

IF  (IPLOT(JM) .GE.5)  GOTO  105 

ITOT=rrOT 

GO  TO  125 

IF  (LOGTYP(.JM)  ) 120,115,110 
ANUX(IK)=XARRAY(IJ) 

ANUY  ( I K)  = (LOG  10 ( YARRAY  ( IH ) ) - DG 10  ( Y.ARRAY  ( 1 TO  1+  1 ) ) ) / 
lYARRAY(IT0T+2) 
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GO  TO  125 

ANUX ( I K ) = ( LOG  1 0 ( XARRAY ( IJ ) ) - LOG 1 0 ( XARR AY ( I T0T+ 1 ) ) ) / 
XARRAY(IT0T+2) 

ANUY (IK) = (LOGIC (YARRAY(IJ))-LOG10(YARRAY(ITOT+1)))/ 
YARRAY(IT0T+2) 

GO  TO  125 

ANUX ( I K)  = (LOG 10 ( XARRAY ( IJ ) ) - LOG 10 ( XARRAY  ( IT0T+ 1) ) ) / 
XARRAY(ITOT+2) 

ANUY(IK)=YARRAY(IJ) 

GO  TO  125 

IK=IK+1 

CONTINUE 

XPTS (IK) =XARRAY (ITOT+ 1) 

XPTS ( I K+ 1 ) =XARRAY ( ITOT+  2 ) 

YPTS (IK) =Y ARRAY (IT0T+ I) 

YPTS (IK+I) =Y ARRAY (ITOT+2) 

IF  (IPL0T(JM).GE.5)  GO  TO  135 
GO  TO  155 

IF  (LOGTYP(JM))  150,145,140 
ANUX ( I K) =XARRAY ( I TOT+ 1 ) 

ANUX ( I K+ 1 ) =XARRAY ( IT0T+ 2 ) 

ANUY(IK)=0.0 
ANUY(IK+1)=1.0 
GO  TO  155 
AlWX(IK)=0. 

ANUX(IK-*-l)=l. 

ANUY(IK)=0. 

ANUY(IK+1)=1. 

GO  TO  155 
ANUX(IK)=0. 

ANUX(IK+1)=1. 

ANUY ( I K) =YARRAY ( IT0T+ 1 ) 

ANUY ( I K+ 1 ) =YARRAY ( ITOT+  2) 

GO  TO  155 
JN=JN+NPTS (JM) 
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TF  (IPL0T(.TM).EQ.3.0R.  IPLOT(JM)  .BJ.O.OR.  IHIST(JM)  .RQ.IUN) 
EO  TO  220 

IF  (IPLOTIJM)  .IiQ.5)  GO  TO  loO 

GO  TO  185 

JZ=3 

ANUHXll)=ANUX(l) 

/VNUHYC1)=ANUY(1) 

ANU11X(2)=ANUX(:) 

ANUI(Y(iJ«ANUV(JJ 
X:=NPTS(JM)-1 
DO  U>5  i:=Z,KZ 
ANUHXCJZ) =ANUX(IZ) 

ANUllY  (JZ)=ANUY(1Z+1) 

JZ=JZ*1 

•ANUHXCJZ)  =ANUX11Z+1) 

ANUllY (JZ)=ANUY  (IZ+1) 

JZ=JZ+1 
CONTI  NUi; 

IF  (LOGTYPCJM))  180,175,170 
ANUI IX  ( JZ)  =XARRAY  { ITOT-t- 1 ) 

ANUIIXCJZ+  1)=XARRAY  (ITOT+Z) 

ANUllY  (JZ)  =0.0 

•ANUI  IY(JZ+1)  = 1.0 

GO  TO  105 

•ANUHXCJZ)  =0.0 

ANUHXCJZ+1)=1.0 

ANUHYCJZ)=0.0 

ANUllY  CJZ+1)  = 1.0 

GO  TO  105 

A.NUHXCJZ)=0.0 

•ANUHX(JZ+1)  = 1.0 

ANUHY  CJZ) =YARRAY  C ITOT+ 1 1 

ANUllY  CJZ+  1)=YARRAY  CITOT+Z) 

GO  TO  105 
J=3 

HIST0XC1)=XPTS(1) 

HIST0YC1)=YPTSC1) 

HISTOXCZ)=XPTSCZ) 

HISTOYCZ)=YPTSCZ) 

KK=NPTS (JM) - 1 
DO  100  I=Z,KK 
HISTOXCJ)=XPTSCI) 

H1ST0YCJ)=YPTSCI+1) 

J=J+1 

HISTOXCJ)=XPT.SCI-*-I) 

H1ST0YCJ)=YI’T.SCI  + 1) 

J=J>1 

CONTINUK 
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HI STOX ( J) tXARRAY ( ITOT +1 ) 

HISTOXCJ* 1) =XARRAY (ITOT^J) 

HISTOY (J) =YARRAY ( ITOT* 1 ) 

HISTOY  (J-^  1)  =YARRAY  (IT0T*2) 

195  NPT=CNPTS(JM)*2)-2 

GO  TO  (200,205,255,210,215,255),  L 

200  CALL  LINE  (HlSTOX.lllSTOY  ,NPT,  1 , LINTYP  (JM) , INTEQ  (JM) ) 

GO  TO  255 

205  CALL  UASHL  (11IST0X,IIIST0Y,NPT,  1) 

GO  TO  255 

210  CALL  LGLIN  (HISTOX, HISTOY, NPT, 1 , LINTYP (JM) , INTEQ (JM) , LOGTYP (JM) 
1 ) 

GO  TO  255 

215  CALL  DASHL  (ANUHX, ANU1IY,NPT,  1) 

GO  TO  255 
C 

C...  SECTION  FOR  PLOTTING  REGULAR  DATA  (NO  HISTOGRAMS). 

C 

C IPLOT=l STRAIGHT  LINE;LINE/\R  AXES 

C IPLOT=2  DASH  LINE: LINEAR  AXES  ONLY 

C IPLOT=3  SMOOTH  LINE: LINEAR  AXES 

C IPL0T=4 STRAIGHT  LINE : LOG- LOG, SEMI - LOG  AXES 

C IPLOT=5  DASH  LINE: LOG-LOG, SEMI-LOG  AXES 

C IPLOT=0 SMOOTH  LINE : LOG- LOG, SEMI - LOG  AXES 

C 

220  GO  TO  (225,230,235,240,245,250),  L 

225  CALL  LINE  (XPTS, YPTS, NETS (JM) , I , LINTYr(JM) , INTEQ (JM) ) 

GO  TO  255 

230  CALL  DASHL  (XPTS,YPTS,NPTS (JM) , 1) 

GO  TO  255 

235  CALL  FLINE  (XPTS  ,YPTS , -NETS  (JM)  , 1 , LINTl'P  (JM) , INTEQ  (JM) ) 

GO  TO  255 

240  CALL  LGLIN  (XPTS, YETS, NPTS (JM) , 1, LINTYP (JM) , INTEQ (JM) , 

1 LOGTYP (JM)) 

GO  TO  255 

245  CALL  DASHL  (ANUX,ANUY,NETS (JM) , 1) 

GO  TO  255 

250  CALL  FLINE  (ANUX,ANUY, -NETS (JM) , 1 , LINTYP (JM) , INTEQ (JM) ) 

255  CONTINUE 

CALL  PLOT  (XAXIS,0.0,3) 

CALL  PLOT  (SKIP, 0.0, -3) 

RETURN 

END 

SUBROUTINE  READR  (NPASS,t4ASK) 

THIS  SUBROUTINE  HANDLES  THE  INTERACTIVE  I/O  - - READS  IN  PLOT 
INSTRUCTIONS  AS  WELL  AS  THE  DATA. 

COMMON  /INPUT/  IXAXIS (5) , lYAXIS (5) , ITITLl (5) , ITITL2 (5) , ITITL3 (5) , 
1ITITL4(5).ITITL5(5),ITYPE.XAXIS,YAXIS,FACT,ITITLE,Y2,X1,IGRU), 
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2YINC,XINC,N01J,N1N,ILINI:S,LU1;V,NS1-TS,111N,11IY,1HB,  lUO 
C 

CONWON  /DATA/  XARRAY  (900) , YARRAY  (900)  ,11ISTOX(900)  .IlISTOY  (900) , 
IXPTS  (900)  ,YPTS(900)  ,ANUX(900)  ,ANIJY  (900)  ,ANU1IX  (900)  .ANIJUY  (900) 

C 

CONMON  /DRAW/  NPTS  (20) , I PLOT (20) , I11IST(20)  ,LINTYP(20) , INTJ;Q(20)  , 
lLa;TYP(20) .SKIP 

...  niH  FOLLOWING  3 CARDS  ARH  FOR  UNIVAC  MACIIINFS. 

DIMl-NSION  IA(2) 

IA(I)=o110ASG,T 
IA(2)=0II  29. 

IF  ((NPASS)  .GT.  1)  GO  TO  12:> 

100  WRITli  (N0U,50C)0) 

RLAl)  (NIN,  5030, F.RR=*  105)  LDFIV 
GO  TO  110 

105  WRITL  (NOU,5020) 

GO  TO  100 

C...  •nili  FOLLOWING  4 CARDS  ARi:  FOR  UNIVAC  MACllINFS. 

C 

110  IF  (LDLV.LQ.O)  GO  TO  115 
GO  TO  120 

115  CALL  LRTRAN  (o.IA) 

LDi;V  = 29 
C 

C...  niL  FOLLOWING  CARD  WOUU)  BH  DSFD  FOR  CDC  MACllINFS. 

C 

Clio  IF  (LDFV.FQ.O)  LDFV=29 
C 

120  CALL  IN  IT 
125  WRITL  (NOU,5070)  NPASS 
WRITL  (NOU,5080) 

RLAD  (NIN,5050)  MASK 
130  WRITL  (NOU,5090) 

RLAD  (NIN,5030,LRR=135)  XAXIS 
CO  lO  140 

135  WRITL  (N0U,5020) 

(•.0  TO  130 

140  IF  (XAXIS. LQ.O)  XAXIS=9 
145  WRITL  (NOU.SlOO) 

RLAD  (NIN,5030,LRR=150)  YAXIS 
GO  TO  155 

150  WRITL  (NOU,5020) 

GO  TO  145 

155  IF  (YAXIS. LQ.O)  YAXIS=7 
U)0  WRITL  (NOU.SllO) 


RiiAD(NIN,5030,ERR=l65)  FACT 
IF  (FACT.EQ.O.)  FACT=1. 

IF  (FACT.GT.O.)  GO  TO  170 
WRITE  (NOU,5120) 

GO  TO  160 

loS  WRITE  (N0U,5020) 

GO  TO  160 
170  CALL  FACTR 

IF  (MASK.NE.IllY)  GO  TO  175 
WRITE  (NOU,5010) 

READ  (NIN,5040)  ITITLl 
GO  TO  235 

175  WRITE  (NOU,5130) 

READ  (NIN,5050)  ITITUi 

IF  (ITITLE.EQ.IHN.  OR.  ITITLE.EQ.IIIB. OR.  ITITLl-;. EQ.IUO)  GO  TO  235 
180  WRITE  (N0U,5140) 

READ  (N1N,5030,ERR=185)  XI, Y2 
GO  TO  190 

185  WRITE  (NOU,S020) 

GO  TO  180 

190  WRITE  (N0U,5150) 

READ  (NIN,5030,ERR=195)  ILINES 
GO  TO  200 

195  WRITE  (NOU,5020) 

1 GO  TO  190 

I 200  WRITE  CNOU,51bO)  ILINES 

'i  DO  230  1 = 1,  ILINES 

GO  TO  (205,210,215,220,225),  I 
, 205  READ  (NIN,5040)  ITITLl 

■ GO  TO  230 

1 210  READ  (N1N,5040)  ITITU 

; GO  TO  230 

215  READ  (NIN,5040)  rriTL3 

GO  TO  230 

220  READ  (NIN,5040)  rriTL4 

i'i  GO  TO  230 

f;  225  READ  (NIN,5040)  ITITL5 

I-  230  CONTINUE 

235  WRITE  (NOU,5170) 

READ  (NIN,5040)  IXAXIS 
WRITE  (NOU,5180) 

READ  (NIN,5040)  lYAXIS 
IF  (NIASK.EQ.  IIIY)  GO  TO  255 
WRITE  (N0U,5190) 

READ  (NIN,5050)  IGRID 

IF  (IGRID.EQ.IHN.OR.  IGRID.  EQ.  IIIB. OR.  IGRID.  IiQ.  1110)  GO  TO  255 
240  WRITE  (NOU,5200) 

READ  (NIN,5030,ERR=245)  LGRID 
GO  TO  250 
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145  WRITI-  CNOU,i>020) 

GO  TO  ’40 

250  XINC=10. 

YINC=10. 

IF  (LGRID.lt. 0)  XINC=0. 

IF  (LGRIU.GT.O)  YINC=0. 

255  WRITL  (NOU,5210) 

RtAU  (NIN,5030,LRR=2o0)  ITYPL 
GO  rO  205 

2oO  WRITL  (NOU,5020) 

GO  TO  255 

205  WRITL  (NOU,5220)  NPASS 

READ  CNIN,5030.LRR=270)  NSLTS 
GO  TO  275 

270  WRITL  (NOU,5020) 

GO  TO  205 

275  11=0 

DO  375  K=l, NSLTS 

IF  (MASK.LQ.IllY.AND.K.LQ.I)  WRITL  (NOIJ,5000) 
IF  (MASK.LQ.IHY)  GO  TO  3o0 

280  WRITL  (NOD, 52 30)  K 

RLAD  (NIN,5050,LRR=285)  LINTYP(K) 

GO  I’O  200 

285  WRITL  (NOD, 5020) 

GO  TO  280 

200  IF  (LIN'l'YP(K)  .GL.O)  GO  TO  205 
IF  (ITYPL.LQ.I)  IPLOT(K)=I 
IF  (ITYPL.GT. I)  I PLOT (K) =4 
IF  (ITYPL.LQ.2)  L0GTYPCK)=1 
IF  (ITYPL.LQ.3)  L0GTYP(K)=-1 
IF  (ITYPL.LQ.4)  LOGTYP(K)=0 
GO  TO  330 

205  IF  (LINTYP(K) .GT.O)  GO  TO  315 

300  WRITL  (N0U,5240)  K 

RLAD  (NIN,5030.LRR=305)  IPLOT(K) 

GO  TO  310 

305  WRITL  (NOD, 5020) 

GO  TO  300 

310  IF  (ITYPL.GT. I)  I PLOT(K) =I PLOT( K) +3 

IF  (ITYPL.LQ.2)  L0GTYP(K)=1 

IF  (ITYPL.LQ.3)  mGTYP(K)  = -I 
IF  (ITYPL.LQ.4)  LCXITYP (K)  =0 
GO  TO  340 

315  WRITL  (NOD, 5250)  K 

RLAD  (NIN,5030,LRR=320)  IPLOT(K) 

GO  TO  325 
WRITL  (NOU,5020) 

GO  TO  315 


320 


325  IF  CIPL0T(K).EQ.2)  IPL0T(K)=3 

IF  (ITYPU.GT. 1)  IPL0T(K)»IPL0T(K)>3 
IF  (mPE.EQ.2)  L0GTYP(K)»1 

IF  (ITYPE.EQ.3)  L0GTYP(K)»-1  * 

IF  (ITYPE.EQ.4)  LOGTYP{K)=0 
330  WRITE  (NOU,52bO) 

READ  (NIN,5030,ERR=335)  INTEQ(K) 

GO  TO  340 

335  WRITE  (NOU,5020) 

GO  TO  330 

340  IF  (IPL0T(K).EQ.3.0R.IPL0T(K).EQ.0)  GO  TO  360 

345  WRITE  (N0U,5270)  K 

READ  CNIN,5050,ERR=350)  IHIST(K) 

GO  TO  355 

350  WRITE  (NOU,5020) 

GO  TO  345 

355  IF  (IHIST(K).EQ.IIIY)  WRITE  (NOU,528a) 

360  WRITE  (NOU,5290)  K 

KFIX=0 
C 

C...  THE  FOLLOWING  LOOP  READS  THE  (X,Y)  DATA. 

C 

DO  365  1=1,100000 
11=11+1 
C 

C...  UN I VAC  FORMATTED  READ. 

C 

READ  (NIN,S030,END=370)  XARRAY (I I ) , YARRAY (I I ) 

C 

C...  CDC  READ  WITH  2-BRANCH  EOF  CHECK. 

C 

C RE  AD (N I N , 5 030 ) XARRAY ( 1 1 ) , Y ARR AY  ( 1 1 ) 

C IF(EOF(NIN))370, 100 

ClOO  IF(ITYPE.EQ.3.AND.XARRAY(II)  .EQ.0.)ianx=l 

C 

IF  (ITYPE.EQ.3.AND.XARRAYCII) .EQ.O.)  KFIX=1 
IF  (ITYPE.EQ.2.AND.YARRAY(II).EQ.O.)  KFIX=1 
IF  C ITY PE . EQ . 4 . AN D. XARRAY ( 1 1 ) . EQ . 0 . OR . I T\ PE . EQ . 4 . AND . YARRAY ( 
1II).EQ.0.)  KFIX=1 
365  CONTINUE 

370  1=1-1 

11=11-1 
NPTS(K)=I 
375  CONTINUE 

IF  (KFIX.EQ. 1)  CALL  FIXUP 
RETURN 
C 

C *****FORMAT  STATEME.NTS***** 

C 

5000  FORMAT  (/51H***  NOTE;  FOR  A 3-0  PLOT  ONLY  A SOLIO-STRAI QIT  LINE,/, 
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14111  IS  AVAILABLE  FOR  ALL  DATA  SHTS.) 

5010  FORMAT  (/43H  PLOT  TITLE  ( ONE  LINE,  25  QIARACTliRS  )?,/, 

125H 25  CHAR  - ) 

5020  FORMAT  (/Soil***  INPUT  ERROR;  RE-ENTER  LAST  LINE.) 

5030  FORMAT  () 

5040  FORMAT  (5Ab) 

5050  FORMAT  (Al) 

5060  FORMAT  (/38H  LOGICAL  UNIT  NUMBER  FOR  PLOT  TAPE?,/, 

12411  DEFAULT  = UNIT  20  .) 

5070  FORMAT  (/38H***  THE  FOLLOWING  PERTAIN  TO  GRAPH  NO.,  12,/) 

5080  FORMAT  (4011  DO  YOU  WANT  A 3-DIMENSIONAL  PLOT  --  YES  OR  NO, 

11711  (DEFAULT  = 2-U)?) 

5000  FORMAT  (/43H  X AXIS  LENGTH  (INDIES)  --  DEFAULT  =0  .) 

5100  FORMAT  (/43H  Y AXIS  LENGTH  (INCHES)  --  DEFAULT  =7  .) 

5110  FORMAT  (/30H  PLOT  SIZE  SCALING  FACTOR  --  DEFAULT, 1411  = FULL  SCA 

ILE.) 

5120  FORMAT  (/43H***  INPlfl’  ERROR:  SCALING  FACTOR  MUST  BE  NON, 
llOH-NEGATIVE, ,/,23H  RE-ENTER  LAST  LINE.) 

5130  FORMAT  (/51H  DO  YOU  WANT  A TITLE  BLOCK  --  YES  OR  NO  (DEFAULT, 
11311  = NO  TITLE)?) 

5140  FORMAT  (/51H  GIVE  X,Y  COORDINATE  (INCHES)  FOR  PLACEMF.NT  OF,/, 
14X,37irniE  UPPER  LEFT  CORNER  OF  TITLE  BLOCK.) 

5150  FORMAT  (/43H  YOU  NOW  HAVE  RCXIM  FOR  5 LINES  OF  PRINT,, 

12411  25  aiARACTERS  PER  LINE.  ,/,4X,27HHOW  MANY  LINES  DO  YOU  NEED?) 

51o0  FORMAT  (/24H  ENTER  SCRIPT  FOR  THE, 12,911  LlNE(S):,/, 

125H 25  CHAR ) 

5170  FORMAT  (/38H  X AXIS  LABEL  (25  CHARACTER  LIMIT)?,/, 

12511  25  QIAR  ) 

5180  FORMAT  (/38H  Y AXIS  LABEL  (25  CHARACTER  LIMIT)?,/, 

125il 25  aiAR ) 

;.iyO  FORMAT  (/441I  DO  YOU  WANT  A GRID  --  YES  OR  NO  (DEFAULT, 

11211  = NO  GRID)?) 


5200  FORMAT  (/4t>ll 
13211  -1=  HORI 

2LKTU;AL  GRID 
FORMAT  (/34H 
IX,  LINEAR  Y,/ 
2Y,/,6X,  16114  a 
FORMAT  (/38H 
FORMAT  (/39il 
14711  0=  POINTS 
24 811 +N a POINTS 
32511  EVERY 


5210 


5220 

5230 


41NTED  AT,/,4X 
:>240  lORmr  (/33H 
119111  a SOLID, 
2U)113  a SOLID, 
5250  FORMAT  (/33H 
119111  a SOLID, 
3260  FORMAT  (/  loH 


S1;LECT  ONi;  of  the  following  for  your  GR1D:,/,4X, 
ZONTAL  GRID  LINES  0NLY,/,4X,4(H1  0=  HORIZONTAL  AND  V 

LINES, /,4X, 3011  +1=  VERTICAL  GRID  LlNliS  ONLY) 

HOW  ARE  Tin:  AXES  TO  BE  SCALED? ,/, OX  ,2211 1 = LINEAR 
,6X,  19112  a LlNli/XR  X,  HX;  Y,/,6X,  19113  = LOG  X,  LINliAR 
LOG  X,  LOG  Y) 

HOW  MANY  SETS  OF  DATA  ON  GRAPH  NO.,  12,  111'’) 
LINE/SYMBOL  COMBINATION  FOR  SET  NO.  , 12,  111; , /,4X, 
CONNliCTED  BY  LINE, NO  SYMBOLS  PR  INTED , /,  4X , 

CONNECTED  BY  LINE,  SYMBOLS  PRINTED  AT,/,4X, 

N-ni  DATA  POINT, /,4X,37H-N=  NO  LINE  DRAWN,  SYMBOLS  PR 
511  EVERY  N-ni  DATA  POINT) 

WIL\  r TYPE.  OF  LINE.  FOR  SF.T  NO.  , 1 2 , 111? , / , 6X  , 

STRAKJIT  ,/,6X,2{MI2  = DASHED,  STRA I GUT  ,/,6X, 

SMOOni) 

WHAT  TYPE  OF  LINE  FOR  SET  NO.  , 1 2 , 111? , / , 6X  , 

STRAIGHT  ,/,6X,  17112  a SOLID, SMOOTH  ) 

INTEGER  EQUIVALENT  OF  CALCOMP  SYMBOL  TO  BE,/,4X, 
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132HPRINTEU  AT  EACH  N-Tll  DATA  POINT?) 

5270  FORMAT  (/39H  DO  YOU  WANT  A HISTOGRAM  FOR  SET  NO. , 12, 111, , /, 4X , 
llOHYES  OR  NO?) 

5280  FORMAT  (/50H***  CAUTION;  FOR  A HISTOGRAM, YOU  NEED  ONLY  SPECIFY,/, 
14X,57H  (1)  A LOWER  AND  UPPER  LIMIT  FOR  THE  FIRST  BIN,  THEN,/, 

24X,54H  (2)  AN  UPPER  LIMIT  ONLY  FOR  EAQl  BIN  'niEREAFTl';R. ) 

5290  FORMAT  (/35H  ADD  YOUR  (X,Y)  DATA  FOR  SET  NO. , 12,  111:  ,/,4X, 

I17HEND  WITH  @EOF  .) 

C 

END 

SUBROUTINE  TITLE 
C 

C THIS  SUBROUTINE  WILL  DRAW  A TITLE  IN  A 3.5  IN.  WIDE  BY  1.5  IN.  HIGH 
C AREA  ANYWHERE  ON  THE  GRAPH  (AS  DEFINED  BY  THE  USER).  THERE  IS 
C SPACE  AVAILABLE  FOR  5 LINES  OF  CHARACTERS,  25  CHARACTERS/LINE. 

C IF  SWITQl  IN  DATA  STATEMENT  IS  SET  TO  1,  SYMBOLS  WILL  BE  PRINTIU) 

C AT  THE  BEGINNING  OF  EACH  TITLE  LINE  CORRESPONDING  TO  niOSE 
C REQUESTED  FOR  PLOTTING  YOUR  DATA  IN  SUBROUTINE  RE.ADR  (INTEl)). 

C 

CONWON  /INPUT/  IXAXIS(5),IYAXIS(5),ITITL1(5),ITITL2(5), 
IITITL3(5),ITITL4(5)  ,ITITL5(5),ITYPE,XAXIS,YAXIS,FACT,rriTLi:,Y2,Xl, 
2IGRID,YINC,XINC,NOU,NIN,ILINES,LDEV,NSETS,IHN,IHY,IHB,  HR) 

C 

COMMON  /DATA/  XARRAY  (900)  ,YARRAY  (900)  ,1IISTOX(900)  ,11IST0Y  (900)  , 
IXPTS  (900)  ,YPTS  (900)  ,ANUX(900)  ,ANUY  (900) , ANU11X(900) , ANUHY  (900) 

C 

COMMON  /DRAW/  NPTS(20)  , IPLOT(2l)) , IHIST(20) , LINTYP(20) , INTEQ(20) , 
1LOGTYP(20),SK1P 
C 

DATA  SWITCH/0./ 

C 

SY=Y2-.25 

SX=XI+.25 

IF  (SWITCH. EQ. 1. ) CALL  SYMBOL  (SX,SY+.05,0. 14, INTEQ(I) ,0.0, -1) 

CALL  SYMBOL  (SX,SY,0. 14,ITITL1,0.0,25) 

IF  (ILINES.lt. 2)  GO  TO  1 
SY-SY-.25 

IF  (SWITCH. EQ. 1. ) CALL  SYMBOL  (SX,SY^ . 05,0 . 14 , INTIU)  (2) ,0 . 0, -1 ) 

CALL  SYMBOL  (SX,SY,0. 14, ITITL2,0.0,25) 

IF  (ILINES.lt. 3)  GO  TO  1 
SY-SY-.25 

IF  (SWITCH.EQ.l.)  CALL  SYMBOL  (SX,SY+ . 05 ,0. 14 , INTEQ (3) , 0. 0. - 1 ) 

CALL  SYMBOL  (SX,SY,0. 14, ITITL3,0 . 0, 25) 

IF  (ILINES.lt. 4)  GO  TO  1 
SY=SY-.25 

IF  (SWITCH.EQ.l.)  CALL  SYMBOL  (SX,SY+. 05 ,0 . 14, INTEQ (4) ,0 .0, -1 ) 

CALL  SYMBOL  (SX,SY.O. 14,rriTL4,0 .0,25) 

IF  (ILINES.lt. 5)  GO  TO  1 
SY=SY-.25 

IF  (SWITCH.EQ.l.)  CALL  SYMBOL  (SX,SY+.  05,0. 14,  INTi;Q(5)  ,0. 0, -1 ) 

CALL  SYMBOL  (SX,SY,0 . 14 , 1TITLS,0.0, 25) 

1 RETURN 

END  ^ 
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I*N*T*E*R*P*0*L 

4 

A general  purpose  2-0IMENSI0NAL/3-DIMENSI0NAL  PLOTTING  ROUTINE 

R.M.  SCHMENK 
BRL  - APG*MO 

THIS  IS  THE  UNIVAC  1108  VERSION-1,  / 27  JAN  78  / 

MAIN  ROUTINE  FOR  * INTERPOL* 

THE  ARRAYS  FOR  THE  (XtY)  DATA  ARE  ARBITRARILY  SET  TO  900 
IN  COMMON  data,  THE  NUMBER  OF  SETS  OF  DATA  PER  GRAPH 
IS  ARBITRARILY  SET  TO  20  VIA  COMMON  DRAW, 

THE  DIMENSION  STATEMENT  (ICHK)  , BELOW  » WAS  ARBITRARILY  SET 
TO  100  — I,E,,  AN  ARBITRARY  NUMBER  OF  GRAPHS  PER  RUN, 

IT  IS  SIMPLY  USED  TO  CHECK  VERTICAL  AXIS  LENGTHS  FOR  PAPER 
SIZE  RESTRICTIONS, 

,,,  THE  NEXT  CARD  IS  FOR  CDC  MACHINES, 

PROGRAM  MAIN( lNPUT,0UTPUT»TAPE5aINPUT.TAPE6=0UTPUT,TAPE29) 

COMMON  /INPUT/  IXAXIS(5) tlYAXlS(5) tlTITLl (5) tlTITLZCS) *ITITL3(5) , 
1ITITLA(5) flTITLSCS) , I TYPE t XAX I S . Y AX  IS .F ACT  * I T I TLE * Y2 . X 1 . IGR ID . 
2YINC»XINCtNOU*NIN*ILINES»LDEV»NSETS»IHN*IHY»lHB.IHO 
C 

COMMON  /DATA/  XARRAY ( 900 ) ♦ YARHAY ( 900 ) ♦ HI STOX ( 900 ) ♦ HI STOY ( 900 ) ♦ 
IXPTS (900) »YPTS (900) ,ANUX (900) ♦ ANUY (900) »ANUHX (900) ,ANUHY (900) 

C 

COMMON  /DRAW/  NPTS ( 20 ) * IPLOT ( 20 ) ♦ I HI  ST (20) *L INTYP ( 20 ) ♦ INTEQ ( 20 ) » 
1LOGTYP(20) tSKIP 
C 

DATA  N0U»NlN/6,5/ 

DATA  IHN,IHY*lHe,lHO/lHN,lHY*lH  »1H0/ 

C 

DIMENSION  ICHK(IOO) 

WRITE  (NOU*5000) 

SKIPsO, 

NPASS=1 

100  CALL  READR  (NPASS*MASK) 

CALL  AXES  (MASK) 

IF  (mask,EO,IHY)  60  TO  105 
IF  (IGRID,EO,lHY)  CALL  NU6RID 

IF  (Ititle,eo,ihy)  call  title 
CALL  PLOTR 
GO  TO  no 

105  call  hide  (YSTEP) 

YAXIS=YAXIS*YSTEP 
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110  ICHK  (NPASS)  *YA)(IS*FACT 
InRITE  (NOU*5010)  NPASS 
READ  (NIN«5020)  IASK 
IF  (lASK.EQ.IHN)  GO  TO  115 
NPASSsNPASS4>1 
GO  TO  100 

115  CALL  PLOT  (SKIP,0.0»999) 

WRITE  (NOU*5030)  NPASS»LDEV 
IPRINT*0 

00  125  laltNPASS 
IF  (ICHK(I)-ll)  125tl25.120 
120  IPRINT»1 

125  CONTINUE 

IF  (IPRINT.EQ.l)  WRITE  (NOU*5050) 

IF  tIPRINT.EQ.O)  WRITE  (NOU*5040) 

WRITE  (NOUf5060) 

*****F0RMAT  STATEMENTS***** 

5000  FORMAT  (/47H  **  INTERPOL.*  UNIVAC  1108  VERSION-1  / 27  JAN  78*/, 
135H  FOLLOW  PRINTED  INSTRUCTIONS,  ) 

5010  FORMAT  (/22H  **  PLOT  FOR  GRAPH  NO.,I2fllH  COMPLETED,*// 

lAlH  >>  DO  YOU  WANT  TO  MAKE  ANOTHER  GRAPH  — tlOHYES  OR  NO<) 
5020  FORMAT  (Al) 

5030  FORMAT  (/3M  ***I2*22H  6RAPH(S)  COMPLETED  — ,/*4X*9HPL0T  FILE* 
118H  WRITTEN  ON  UNIT  ••*I2*2H'',) 

5040  FORMAT  {/39H  *•  PLOT(S)  WILL  FIT  ON  "NARROW**  PAPER,) 

5050  FORMAT  (/43H  **  PLOT(S)  WILL  NOT  FIT  ON  "NARROW"  PAPER. ,/*4X* 
116HUSE  "WIDE"  SIZE,) 

5060  FORMAT  (/20H  **  END  OF  RUN /) 

END 

SUBROUTINE  LOOKUP  (X*XTBL*J) 


THIS  SUBROUTINE  IS  CALLED  FROM  HIDE  AS  AN  EFFICIENT  METHOD 
OF  table  LOOKUP  (STOLEN  FROM  TRACOR  3-D  PLOT  PACKAGE), 

DIMENSION  XTBL(l) 

Ja2 

1 IF  (XTRL(J)-X)  2*3*4 

2 JaJ*l 
GO  TO  1 

3 RETURN 

4 J=J-1 
RETURN 
END 

SUBROUTINE  FIXUP 

THIS  routine  sets  ANY  DATA  WHICH  IS  TO  BE  LOGARITHMICALLY  SCALED 


1^ 


a. 


i 
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C AND  WAS  FOUND  TO  BE  EQUAL  TO  ZERO  TO  THE  MINIMUM  OF  THE  RESPECTIVE 
C X OR  Y ARRAY,  THIS  PREVENTS  AN  ABORTED  RUN.  THE  USER  IS  TOLD. 

C 

COMMON  /input/  IXAXIS(5) »IYAXIS(5) ,ITITLl (5) » ITITL2 (S) ♦ ITITL3 «5) , 
IITITLA (5) .ITITLSCS) »ITYPE»XAXlS*YAXlS*FACT,ITITLE.Y2»Xl*IGRIDf 
2YINC«XINC*NOU*NlN*ILINEStLDEV*NSETS*IHN*lHY*lHB.lHO 
C 

COMMON  /DATA/  XARRAY (900) t YARRAY (900) tHISTOX (900) ^HISTOY (900)  ♦ 
1XPTS(900) »YPTS(900) «ANUX (900) ♦ANUY (900) .ANUNX (900) ,ANUHY (900) 

C 

COMMON  /DRAW/  NPTS(20) *IPLOT(20) tlHlST(20) .LINTYP ( 20 ) * INTEQ ( 20 ) . 
ILOGTYP(20) .SKIP 
C 

DO  1 lal.NSETS 
ITOT*ITOT*NPTS(I) 

1 continue 

XMIN=1 .0E*29 
YMINsl .0E^29 

DO  2 I=l,ITOT 
XMINsAMINl (XMIN.XARRAY(  I) ) 

YMINsAMINI (YMIN.Y ARRAY ( I ) ) 

2 CONTINUE 

IF  (XMIN.EQ.O.)  XMIN=XMIN^.0001 

IF  (YMIN.EQ.O.)  YMINsYMlN^.OOOl 

L=ITYPE-l 
GO  TO  (3.5.7)  . L 

3 DO  4 lal.ITOT 

IF  (YARRAY(I)  ,EQ.O.)  YARRAY ( I ) »YMIN 

4 CONTINUE 

GO  TO  9 

5 no  6 Isl.ITOT 

IF  (XARRAY(I)  .EQ.O.)  XARRAY ( I ) sXMlN 

6 CONTINUE 

GO  TO  9 

7 DO  8 I=l,ITOT 

IF  (XAPRAY(I)  ,EQ.O.)  X ARRAY ( I ) *XM IN 
IF  (YARRAY(I)  .EQ.O.)  YARRAY ( I ) sYMiN 

8 CONTINUE 

9 CONTINUE 
WRITE  (NOU.IO) 

RETURN 

format  STATEMENTS***** 

10  FORMAT  (/51H  **  INPUT  WARNING:  ZER0»S  ENCOUNTERED  ON  LOGARITHMI. 
119HCALLY  SCALED  DATA— ./.36H  VALUES  WERE  RE-SET  TO  MINIMUM  OF, 
233H  RESPECTIVE  ARRAY  (XMIN  OR  YMlN).) 

C 

END 


30 


oooo  o o o n ooooo 


SUBROUTINE  INIT 

THIS  SUBROUTINE  INITIALIZES  THE  PLOTTER.  THE  (X»Y)  COORDINATES 
IN  THE  DATA  STATEMENT  SET  THE  POSITION  OF  THE  6RAPH»S  ORIGIN 
ON  THE  PAPER  (I.E.»  LOWER  LEFT  HAND  CORNER). 

COMMON  /input/  IXAXIS(5) .IYAXIS(5) .ITITLl (5) ♦ITITL2(5) .ITITLGCS) . 
1ITITLA(5) .ITITL5(5) f ITYPE.XAXISfYAXIS.FACT,ITITLE»Y2»Xl.IGRID. 

2yinc*xinc»nou*nin.ilines»ldev.nsets«ihn»ihy» ihb.iho 

COMMON  /DATA/  XARRAY (900 ) f YARRAY ( 900 ) .HI STOX ( 900 ) .HI STOY (900 ) » 
1XPTS(900) .YPTS(900) .ANUX (900) .ANUY (900 ) . ANUHX ( 900 ) .ANUHY (900) 

COMMON  /DRAW/  NPTS (20) . IPLOT (20) . IHIST (20) .LINTYP(20) ♦INTEQ(20) » 
1L06TYP(20) .SKIP 

DATA  X,Y/3. 0.2.0/ 

CALL  PLOTS  (I0UFF.NLOC.LDEV) 

CALL  PLOT  (0.0. -36. 0.-3) 

CALL  PLOT  (X.Y.-3) 

RETURN 

END 

SUBROUTINE  FACTR 

THIS  SUBROUTINE  CHANGES  THE  LENGTH  OF  ALL  PEN  MOVEMENTS  BY  A 
USER  SUPPLIED  SCALING  FACTOR  (SU0R.  READR).  1.0  IS  FULL  SCALE. 

COMMON  /input/  IXAXIS(5)  .IYAXIS(5)  .ITITLKS)  ♦ITITL2(F.)  . 

IITITL3(5) .ITITL4(5) . IT ITL5 ( 5 ) . I TYPE . XAXI S . YAX I S.FACT . I T I TLE . Y2 . X I . 
2IGRID«YlNC.XlNC.N0U.NlN»ILINES.L0EV.NSETS.IHN.IHy .IHB.IHO 

COMMON  /DATA/  XARRA Y ( 900 ) . YARRA Y ( 900 ) . HISTOX ( 900 ) . HISTO Y ( 900 ) . 
IXPTS (900) .YPTS (900) .ANUX (900) .ANUY (900) .ANUHX (900) .ANUHY (900) 

COMMON  /DRAW/  NPTS (20) .IPLOT (20) . IHIST (20 ) .LINTYP (20) .INTEQ(20) . 
1LOGTYP(20) .SKIP 

CALL  factor  (FACT) 

RETURN 

END 

SUBROUTINE  NUGRID 

THIS  SUBROUTINE  DRAWS  GRID  LINES  AS  REQUESTED  BY  THE  USER.  IT 
WILL  NOT  DRAW  THROUGH  THE  TITLE.  ZINC  IS  THE  INCREMENT  BY  WHICH 
THE  PLOT  PEN  CHECKS  TO  SEE  IF  IT  IS  WITHIN  THE  TITLE  BOUNDARIES 
SET  UP  BY  SUBROUTINE  TITLE. 

COMMON  /input/  IXAXIS(5)  .IYAXIS(5)  .ITITLKS)  .ITITL2(5)  .ITITL3(S)  . 


non 


IITITL4(5) *ITITL5(5) » ITYPE»XAXIS»YAXIS»FACT, ITITLE*Y2*X1 * IGRIUt 
2YINC»XINC»NOU*NIN*ILINES*LDEV»NSETS*IHN*IHY,IH8»IHO 

COMMON  /DATA/  X ARRAY (900 ). YARRAY (900 ) t HISTOX (ROO ) *HISTOY (900 ) • 
IXPTS (900) .YPTS (900) .ANUX (900) »ANUY (900) ^ANUNX (900) ,ANUHY(900) 

COMMON  /DRAW/  NPTS(20) ♦IPLOT(20) •IHIST(20) ♦LINTYP(20) tINTEQ(20) ♦ 
IL06TYP(20) fSKIP 

DATA  ZINC/0.1/ 

IF  (ILINES.EQ.O)  60  TO  125 

GO  TO  (100*105.110.115*120).  ILINES 

YlaY2-.3 

60  TO  125 

YlaY2-.6 

60  TO  125 

Yl=Y2-,9 

60  TO  125 

Y1=Y2-1.2 

60  TO  125 

Y1»Y2-1.5 

CONTINUE 

X2sX1.3,75 

DRAWS  VERTICAL  GRID  LINES. 

ITOT=0 

DO  130  I»1.NSETS 
ITOTsITOT.NPTS( I) 

CONTINUE 

IF  (XINC.EQ.0.0)  60  TO  160 
IF  (ITYPE.EQ.1.0R.ITYPE.EQ.2)  XINC=1. 

IF  (ITYPE.EQ.3.0R.ITYPE.EQ.4)  XINC*! ./XARRAY ( ITOT.2) 

X=XINC 

ICOUNT*XAXIS/XINC 
JCOUNTsYAXIS/ZINC 
DO  155  lal.ICOUNT 
YaO.O 

CALL  PLOT  (X.Y.3) 

DO  150  Jal.JCOUNT 
YsY.ZINC 

IF  (X.6E.XI.AN0.X.LE.X2)  GO  TO  135 
60  TO  lAO 

IF  (Y.6E.Y1.AN0.Y.LE.Y2)  GO  TO  U5 
CALL  PLOT  (X.Y.2) 

GO  TO  150 

CALL  PLOT  (X.Y.3) 

CONTINUE 
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...  DRAMS  HORIZONTAL  GRID  LINES. 


IF  (YINC.EQ.0.0)  60  TO  190 
160  IF  (ITYPE.EQ.1.0R.ITYPE.EQ.3)  YINC=1, 

IF  (ITYPE.EU.2.0R.ITYPE.EQ.A)  YINC=1 ./YARRAY ( I T0T*2) 

YsYINC 

> kcount«yaxis/yinc 

LCOUNTsxAXiS/ZINC 
DO  185  Ksl.KCOUNT 
XsO.O 

CALL  PLOT  (X*Y,3) 

DO  180  Lal*LC0UNT 
XsX^ZINC 

IF  <Y.GE.Y1.AND.Y.LE.Y2)  GO  TO  165 
GO  TO  170 

165  IF  (X.6E.X1.AN0.X.LF.X2)  GO  TO  175 

170  CALL  PLOT  (X.Y.2) 

GO  TO  180 

175  CALL  PLOT  (X,Y*3) 

180  CONTINUE 

Y=Y*yIMC 
185  continue 

190  RETURN 
END 

SUBROUTINE  TITLE 
C 

C THIS  SUBROUTINE  WILL  DRAW  A TITLE  IN  A 3.5  IN.  WIDE  BY  1.5  IN.  HIGH 
C AREA  ANYWHERE  ON  THE  GRAPH  (AS  DEFINED  BY  THE  USER).  THERE  IS 
C SPACE  available  FOR  5 LINES  OF  CHARACTERS*  25  CHARACTERS/LINE. 

C IF  SWITCH  IN  DATA  STATEMENT  IS  SET  TO  1,  SYMBOLS  WILL  BE  PRINTED 
C AT  THE  BEGINNING  OF  EACH  TITLE  LINE  CORRESPONDING  TO  THOSE 
C REQUESTED  FOR  PLOTTING  YOUR  DATA  IN  SUBROUTINE  READP  (INTEQ). 

C 

COMMON  /input/  IX AX  I S ( 5) . I Y AXIS (5) . IT ITL 1 ( 5) . I TI TL2 ( 5) , 

1ITITL3(5) tlTITLACS) . ITITL5 (5) * I TYPE ♦ X AX  I S . YAX I S . F ACT  * I TI TLE » Y2  ♦ X 1 , 
2IGRI0«YINC*XINC*N0U*NIN. ILINES.LDEV.NSETS* IHN. IHY. IHB. IHO 
C 

COMMON  /DATA/  XARRAY ( 900 ) ♦ YARRA Y ( 900 ) . HI STOX ( 900 ) tH I STOY ( 900 ) * 
IXPTS (900) »YPTS (900) « ANUX (900) .ANUY (900) . ANUHX (900) . ANUHY (900) 

C 

COMMON  /DRAW/  NPTS(20) ♦IPLOT(20) *IHIST(20) .LlNTYP(20) tINTEQ(20) t 
1LOGTYP(20) tSKIP 
C 

DATA  SWITCH/0./ 

C 

SYSY2-.25 
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SX=X1>,25 

IF  (SWITCH.EU.l. ) CALL  SYMBOL  (SX ♦ SY* . 05 « 0 . 14 , INTEO ( 1 ) . 0 . 0 ♦-!  ) 
CALL  SYMBOL  ( S X * SY * 0 . 1 4 , ITITL 1 1 0 . 0 » 25) 

IF  (ILINES.LT. 2)  GO  TO  1 
SY=SY-.25 

IF  (SWITCH.EQ.l.)  CALL  SYMBOL  ( SX . SY* . 05 « 0 . 14 , INTEQ ( 2) . 0 . 0 .-1 ) 
CALL  SYMBOL  ( SX , SY » 0 . 1 4 , I T ITL2 t 0 . 0 t 25) 

IF  (ILINES.LT. 3)  60  TO  1 
SYaSY-.25 

IF  (SWITCH.EQ.l.)  CALL  SYMBOL  ( SX » SY* . 05 » 0 . 1 4 , INTEQ ( 3 ) . 0 . 0 . -1 ) 
CALL  SYMBOL  ( SX »SY , 0 . 14 , IT ITL3 t 0 • 0 ♦ 25 ) 

IF  (ILINES.LT. 4)  GO  TO  1 
SYaSY-.25 

IF  (SWITCH.EQ.l.)  CALL  SYMBOL  {SX«SY* . 05 ♦ 0 . 1 4, INTEQ ( 4) . 0 . 0 , -1 ) 
CALL  SYMBOL  ( SX . SY . 0 . 14 . IT ITL4 . 0 . 0 ♦ 25 ) 

IF  (ILINES.LT. 5)  GO  TO  1 
5Y=SY-,25 

IF  (SWITCH.EQ.l.)  CALL  SYMBOL  ( SX ♦ SY* . 05 » 0 . 1 4 , InTEQ ( 5 ) . 0 . 0 ♦ -1 ) 
CALL  SYMBOL  (SX.SY,0.14.ITITL5.0.0.25) 

1 RETDMN 
END 

SUBROUTINE  AXES  (MASK) 

THIS  SUBROUTINE  SCALES  ALL  DATA  AND  DRAWS  THE  APPROPRIATE  AXES 
WITH  LAPELS  AS  WELL  AS  THE  GRAPH'S  BORDER. 

COMMON  /INPUT/  IX AX  IS (5) .IYAXIS(5) flTITLl (5) .ITITL2(5) .ITITL3(5) , 
1ITITL4(5) ♦ITITL5(5) . ITYPE» XAXIS* YAXlSt FACT. ITITLE* Y2» X 1 ♦ IGRID* 
2YINC.XINC.NOU.NIN. ILINES.LDEV.NSETS. IHN.IHY, IH8» IHO 

COMMON  /DATA/  X ARRA Y ( 90 0 ) . Y ARRA Y ( 90 0 ) . H I STOX ( 90 0 ) . H I STO Y ( 90 0 ) . 
1XPTS(900) .YPTS(900) .ANUX(900) .ANUY(900) .ANUHX(900) .ANUHY(900) 

COMMON  /DRAW/  NPTS(20) . I PLOT (20) ♦IHIST(20) ♦LlNTYP(20) .InTEQ(20) . 
1LOGTYP(20) .SKIP 

ITOT=0 

DO  100  I=1.NSETS 
ITOT=ITOT*NPTS(I) 
mo  CONTINUE 

60  TO  (105.110.115.120).  ITYPE 

...  ITYPEsl  — — > LINEAR  X. LINEAR  Y SCALED  AXES. 

105  CALL  scale  ( XARHAY. XAX IS. ITOT. 1 ) 

CALL  scale  (YARRAY.YAXIS.ITOT.1) 

IF  (MASK.EQ.IHY)  GO  TO  125 

CALL  AXIS  (O.O.O.O.IXAXIS.-25.XAXIS.O.O.XARRAY(ITOT*1) . 

IXARRAY ( IT0T*2) ) 
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C 

c • • • 

c 

110 


c 

c ■ • • 

c 

115 


C 

c • • • 

c 

120 


CALL  AXIS  (0.0«0.0tIYAXISf25*YAXlS«90.0tYARRAY(ITOT«l) t 
lYARRAY (IT0T^2) ) 

CALL  PLOT  (0.0,YAXIS.3) 

CALL  PLOT  (XAXIS»YAXIS*2) 

CALL  PLOT  (XAXIS*0.0f2> 

GO  TO  125 

ITYPE»2  — — > LINEAR  X*LO-  Y SCALED  AXES. 

CALL  SCALE  (XARRAY»XAXIS.IT0T.1) 

CALL  SCAL6  ( YARRAY. YAXIS* ITOT» 1 ) 

IF  (MASK.EO.IHY)  60  TO  125 

call  AXIS  (O.OtO.O«lXAXlS«-25f XAXiSfO.OtXARRAY (ITOT«l ) * 
1XARRAY(IT0T*2) ) 

CALL  LGAXS  ( 0 . 0 * 0 . 0 ♦ I YAXI S ♦ 25 1 YAXIS . 90 . 0 » YARRAY ( ITOT* 1 ) » 
1YARRAY(IT0T^2) ) 

CALL  PLOT  (0.n,YAXIS«3) 

CALL  PLOT  (XAXIS*YAXIS»2) 

CALL  PLOT  (XAXISt0,0,2) 

GO  TO  125 

ITYPE=3  LOG  X. LINEAR  Y SCALED  AXES. 

CALL  SCALG  < X ARRA Y » X AX  I S . I TOT ♦ I ) 

CALL  scale  (YARRAY.YAXIS* IT0T» 1 ) 

IF  (MASK.EQ.IHY)  GO  TO  125 

CALL  LGAXS  ( 0 . 0 ♦ 0 . 0 . I X AX  IS ♦ -25 1 XAX I S • 0 . 0 ♦ XARRAY ( ITOT* 1 ) ♦ 
1XARRAY(IT0T>2) ) 

CALL  AXIS  (0.0.0.0«IYAXIS*25#YAXIS.90.0»YARRAY(ITOT*1) ♦ 
1YARRAY(IT0T>2) ) 

CALL  PLOT  (0.0.YAXIS*3) 

CALL  PLOT  (XAXIStYAXIS*2) 

CALL  PLOT  (XAXlS*0.0f2) 

GO  TO  125 

ITYPE=4  — — > LOG  X.LOG  Y SCALED  AXES. 

CALL  SCALG  ( XARRAY ♦ X AX  IS ♦ I TOT ♦ 1 ) 

CALL  SCALG  ( Y ARRA Y . YAX IS ♦ ITOT t 1 ) 

IF  (MASK.EQ.IHY)  GO  TO  125 

CALL  LGAXS  ( 0 . 0 1 0 . 0 ♦ I XAX I S » -25 ♦ XAX I S . 0 . 0 « XARRA Y ( I TOT* 1 ) ♦ 
1XARRAY(IT0T*2) ) 

CALL  LGAXS  ( 0 . 0 1 0 . 0 ♦ I YAX I S * 25 » YAX I S ♦ 90 . 0 ♦ YARRA Y ( I TOT* 1 ) » 
lYARRAY ( IT0T*2) ) 

CALL  PLOT  (0.0.YAXIS*3) 

CALL  PLOT  (XAXIS.YAXIS.2) 

CALL  PLOT  (XAXIS»0.0»2) 

RETURN 

END 


i 
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SUBROUTINE  PLOTR 


THIS  subroutine  SETS  UP  THE  DATA  ARRAYS  WITH  APPROPRIATE  SCALING 
FACTORS  (COMPUTED  IN  SUBR.  AXES) *GENERATES  HISTOGRAMS  IF  REQUESTED 
AND  THEN  PLOTS  ALL  DATA  AS  INPUT  IN  SUBROUTINE  REAOR. 

THIS  ROUTINE  IS  ALSO  NEXT  TO  IMPOSSIBLE  TO  EXPLAIN  OR  FOLLOW! 

COMMON  /INPUT/  IXAXIS(5) « I VAX  IS (5) tlTITLl (5) .ITITL2(5) tlTITLBCS) 
IITITLA(5) »ITITL5(5) ♦ ITYPE f X AX ISt Y AX  I S tF ACT , I T I TLE » Y2* X I » I GRID » 

2yinc*xinc*nou*nin*ilines»loev,nsets»ihn*ihy.ihb*iho 

COMMON  /DATA/  XARRAY (900 ) ♦ YARRAY (900 ) tHISTOX (900 ) ♦HISTOY (900 ) ♦ 
1XPTS(900) *YPTS(900) tANUX (900) »ANUY (900) fANUHX (900) ,ANUHY (900) 

COMMON  /DRAW/  NPTS(20)  t IPLOT(20)  ♦IHIST(20)  *LINTYP(20)  t INTE(J(20)  « 
1LOGTYP{20) fSKIP 

SKIP=XAXIS^4.0 
JN=l 
JO  = 0 
ITOT=0 

DO  100  Isl«NSETS 
ITOT=rTOT«-NPTS(I) 

100  CONTINUE 

C 

DO  255  JM=l*NSETS 
L=IPL0T(JM) 

J0=J0*NPTS ( JM) 

IK=1 

DO  130  IJ=JN,J0 
XPTS(IK)=XARRAY(IJ) 

YPTS(IK)aYARRAY(IJ) 

IF  ( IPLOT ( JM) .GE.5)  GO  TO  105 

ITOT=ITOT 

GO  TO  125 

105  IF  (LOGTYP(JM))  120tll5*110 

110  ANUX(IK)s:XARPAY(IJ) 

ANUY ( IK) = (LOGIC (YARRAY ( I J) ) -LOGIO (Y ARRAY ( ITOT*l ) ) ) / 

1 YARRAY(IT0T*2) 

GO  TO  125 

115  ANUX (IK)s(LOGlO (XARRAY (IJ) ) -LOG  1 0 ( XARRAY ( I TOT ♦ 1 ) > ) / 

1 XARRAY(IT0T42) 

ANUY ( IK) a (LOGIO (YARRAY ( I J) ) -LOG  1 0 ( Y ARRAY ( I TOT> I ) ) )/ 

1 YARHAY(IT0T42) 

GO  TO  125 

120  ANUX ( IK) s (LOGIO (XARRAY ( I J) ) -LOGIC (XARRAY ( ITOT^l ) ) ) / 

1 XARRAY ( IT0T*2) 

ANUY ( IK) sYARRAY ( I J) 


non 


125 

130 


t 

! 

» 

i 

j 

135 

140 


145 


i 

i 

\ 150 

i 


( 


155 


60  TO  125 

IK=IK^1 

CONTINUE 

XPTS(IK)=XARRAY(IT0T*1) 

XPTS(IK*1)=XARRAY(IT0T^2) 

YPTS(IK)=YARRAY{IT0T*1) 

YPTS (IK*!) sYARRAY ( I TOT *2) 

IF  (IPLOT(JM) .GE.5)  GO  TO  135 
60  TO  155 

IF  (LOGTYPCJM))  150»145»140 
ANUX ( IK) sXARRAY ( IT0T*1 ) 

ANUX (IK*1)=XARRAY (IT0T*2) 

ANUY(IK)=0.0 

ANUY( IK*l)sl.O 

60  TO  155 

ANUX(IK)=0. 

ANUX(IK*1)=1. 

ANUY(IK)=0. 

ANUY (IK*1)=1. 

GO  TO  155 
ANUX ( IK) =0. 

ANUX ( IK*l ) =1, 

AMUY (IK)=YARRAY (IT0T*1  ) 

ANUY { IK*l ) =Y ARRAY ( IT0T*2) 

60  TO  155 
JN=JN*NPTS( JM) 


...  SECTION  TO  GENERATE/PLOT  NEw  ARRAYS  FOR  HISTOGRAM, 

IF  (IPLOT(JM) .EQ.3,0R.IPL0T(JM) ,£Q.6.0R, IHIST ( JM) ,£Q, IHN) 
1 GO  TO  220 

IF  ( IPLOT ( JM) ,EQ.5)  GO  TO  160 
60  TO  185 
160  JZ=3 

ANUHX ( 1 ) sANUX  ( 1 ) 

ANUHY(1)=ANUY(1) 

ANUHX (2) =ANUX (2) 

ANUHY (2) =ANUY (2) 

KZ=NPTS(JM)-1 

DO  165  IZ=2»KZ 
ANUHX (JZ) =ANUX (IZ) 

ANUHY (JZ) =ANUY  ( IZ*1  ) 

JZaJZ*l 

ANUHX (JZ) sANUX ( IZ*1  ) 

ANUHY(JZ)sANUY(IZ*l) 

JZsJZ*! 

165  CONTINUE 

IF  {L06TYP(JM))  16fl»175»lT0 
170  ANUHX (JZ) sXARRAY ( IT0T*1 ) 
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ANUHX  (JZ+U  ax  ARRAY  (IT0T*2) 

ANUHY  < JZ) aO.O 
ANUHY(JZ^l)sl.O 
60  TO  195 

175  ANUHX (JZ)aO.O 

ANUHX(JZ*l)al,0 
ANUHY(JZ>=0.0 
ANUHY(JZ^l)al.O 
60  TO  195 

180  ANUHX ( JZ) =0,0 

ANUHX(JZ^l)al.O 

ANUHY (JZ) =Y ARRAY ( ITOT^l) 

ANUHY ( JZ+ I ) aVARRAY ( ITOT^Z) 

60  TO  195 
185  Ja3 

HIST0X(1)=XPTS(1) 

HIST0Y(1)=YPTS(1) 

HISTOX (2)=XPTS(2) 

HISTnY(2)=YPTS(2) 

KusMPTS ( JH) -1 
on  190  I=2*kk 
HrSTOX(J)aXPTS(I) 

HlSTOY(J)aYPTS(I*I) 

J=J*1 

HISTOX{J)aXPTS(I*l) 

HIST0Y(J)=YPTS(I*1) 

JaJ^l 

190  CONTINUE 

HISTOX (J)=XARRAY(IT0T>1) 

HISTOX (J*l)aXARPAY(IT0T42) 

HISTOY(J)aYARRAY(ITOT^l) 

HIST0Y{J*1) ay array (ITOT+2) 

195  NPTa(NPTS(JM)*2)-2 

60  TO  (200»205*255f2l0»215»255) , L 
200  CALL  LINE  (H I STOX . H I STO Y ♦NPT ♦ 1 »L INTYP ( JM ) * INTEQ ( JM ) ) 


60  TO  265 

205  CALL  DASHL  (HISTOX ♦HISTOY*NPT» 1 ) 

GO  TO  255 

210  CALL  LGLIN  (HIST0X»HIST0Y*NPT» 1 ^LINTYP ( JM) , InTEQ ( JM) ,L0GTYP ( JM) 

1 ) 

GO  TO  255 

215  CALL  DASHL  ( ANUHX ♦ ANUHY tNPT  ♦ 1 ) 

GO  TO  255 

...  SECTION  FOR  PLOTTING  REGULAR  DATA  (NO  HISTOGRAMS). 

IPLOTal  -— > STRAIGHT  LINE.'LINEAR  AXES 
IPLOTa?  — — > DASH  LiNEtLINEAR  AXES  ONLY 
IPLOTaS  — — > SMOOTH  LINE '.LINEAR  AXES 


f. 
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IPL0TS4  - — ->  STRAIGHT  L INE ; LOG-LOG ♦ SEMI-LOG  AXES 
IPLOTaS  — — > DASH  LINE ;L0G-L0G»SEMI-L06  AXES 
IPL0T=6  — — > SMOOTH  L INE : LOG-LObt SEMI -LOG  AXES 

220  60  TO  (22S«230*235*240«245*250) * L 

225  CALL  LINE  ( XPTS ♦ YPTS t NPTS ( JM) ♦ I tL INTYP ( JM) ♦ InTEQ ( JM ) ) 

GO  TO  255 

230  CALL  OASHL  ( XPTS * YPTS «NPTS ( JM ) t 1 ) 

60  TO  255 

235  CALL  FLINE  (XPTS*YPTS*-NPTS ( JH) , 1 ♦LINTYP ( JM) ♦ INTEQ < JM) ) 

60  TO  255 

240  CALL  LGLIN  (XPTS* YPTS*NPTS ( JH) • 1 .LINTYP ( JM) ♦ INTEQ ( JM) ♦ 

1 L06TYP(JM)) 

GO  TO  255 

245  CALL  OASHL  ( ANUX ♦ ANUY ♦ NPTS ( JM ) » 1 ) 

GO  TO  255 

250  CALL  FLINE  ( ANUX . ANUY t-NPTS ( JM ) , 1 ,L INT YP ( JM ) ♦ INTEQ ( JM) ) 

255  CONTINUE 

CALL  PLOT  (XAXIS»0.0»3) 

CALL  PLOT  (SKIP*0.0»-3) 

RETURN 

END 

SUaROUTiNE  HIDE  (YSTEP) 

THIS  ROUTINE  PLOTS  THE  3-DIMENSIONAL  GRAPH. 

COMMON  /INPUT/  IXAXIS(5) .IYAXlS(5) . ITITLl (5) »ITITL2(5) .ITITL3(5) 
1ITITL4(5) .ITITL5{5) . I TYPE . XAXI S* YAXIS * FACT  * I T I TLE ♦ Y2 • X 1 . IGRID . 
2YINC.XINC»N0U*NINt ILINES»LDEV»NSETS»IHN»ImY» IHB«IH0 
C 

COMMON  /DATA/  XARRAY ( 900 ) . YARRAY ( 900 ) » H ISTOX ( 900 ) « HI STOY ( 900 ) ♦ 
1XPTS(900) ♦YPTS{900) ♦ANUX(900) »ANUY (900) ♦ANUHX(900) ♦ANUHY(900) 

C 

COMMON  /DRAW/  NPTS ( 20 ) . I PLOT ( 20 ) »IHIST(20)  .LINTYP(20)  ♦INTEQ(20)  . 
1L06TYP(20) .SKIP 
C 

DIMENSION  X (900) «Y (900) .XG (900) .6(900) .XH (900) .H (900) 
equivalence  (KI.IWHICH).  (K2.SL0Pt).  (FNSMl.Zl).  (I66P1.K1),  (M 
1N2) 

DATA  EPSl.MAXOIM.XSTART.YSTART/l.E-9.900.4, ,3./ 
F(XX.XI.YI.XIP1.YIP1)=YI* (XX-XI)*{YIP1-YI)/(XIPI-XI) 

WRITE  (NOU.5000) 

READ  (NINt5030)  IDASH 

YSTEPsYSTART 

SKIP=XAXIS*4.0 

JNsl 

JOsO 

NGaO 

ITOTaO 


100 

C 


105 

1 


110 

1 


115 

1 

1 


120 

125 

130 


135 


140 


DO  100  Isl*NSETS 

ITOT=ITOT*NPTS(I) 

CONTINUE 


00  355  JMsI.NSETS 
NlaNPTS(JM) 

JOajO't'Nl 

IKsl 

DO  130  IJaJNtjO 
ITOTalTOT 

GO  TO  (120*105»110.115)  . ITYPE 

Y(IK)a(LOG10(YAPRAY<IJ) ) -LOGIO ( Y ARRAY ( ITOT* 1 ) ) ) /YARRAY ( I TOT* 
2) 

X(IK)3XARRAY(IJ) 

GO  TO  125 

X(IKJ*(L0610(XARRAY<IJ) )-L06lO(XARRAY(ITOT*l) ) )/XARRAY(ITOT* 
2) 

Y(IK)aYARRAY(IJ) 

GO  TO  125 

X (IK)  = (LOGIC  (XARPAYdJj  ) -LOGl  0 { XARPAY  ( I TOT*  1 ) ) ) /XARRAY  ( I TOT* 
2) 

Y( IK)a(LOG10(YARPAY{IJ) )-LUGlO(YARPAY(ITOT*l) ) ) /YARPAY ( I TOT* 
2) 

GO  TO  125 
X(IK)aXARRAY(IJ) 

Y ( IK) aYARRAY ( IJ) 

IKaIK*l 

CONTINUE 

JNaJN*Nl 

00  135  laJ.Nl 

IF  (X(I-l) .LT.X(I) ) GO  TO  135 
WRITE  <NOU*5010)  X(I-1)*X(I) 

RETURN 

CONTINUE 

IF  (JM.GT.l)  60  TO  205 
NFMSaNSETS 


XNINaXARRAYdTOT*!) 

YMINaYARRAYdTOT*!) 

OELTAXaX ARRAY  dT0T*2) 

DELTAYaYARRAY  dT0T*2) 

IF  (N1*4.LE.MAXDIM)  GO  TO  140 

GO  TO  360 

SIGNal, 

IF  (NG.LT.-I)  SIGNa-1. 

IF  (N6.EQ.-1.0R,NG.EQ.-3)  GO  TO  1a5 
CALL  PLOT  (0.fYSTART*YAXIS*3) 

CALL  OASHP  (XAXIS»YSTART*YAXIS»,03) 
call  OASHP  (XAXIS«YSTAPT» .03) 

CALL  OASHP  (0..YSTART..03) 


40 


CALL  PLOT  (XAXlStYSTART*3) 

CALL  OASHP  (XAXIS*XSTART*0.«.03) 

1A5  CALL  SYMBOL  (2. 0 1 YAX IS^YSTART* . 10  * . Ut I T ITLl * 0 . ,25) 

60  TO  (150tl55.l60.165).  ITYPE 

150  CALL  AXIS  (XSTART.O..IXAXIS.-25.XAXIS.O..XMIN.OELTAX) 

CALL  PLOT  (XSTAPT.0..3) 

CALL  PLOT  <0..YSTART.2) 

CALL  AXIS  (0..YSTART.IYAXIS.25.YAXIS.90..YMIN.DELTAY) 
60  TO  170 

155  call  axis  (XSTART.O..IXAXIS.-25.XAXIS.O..XMIN.OELTAX) 

CALL  PLOT  {XSTART.0..3) 

CALL  PLOT  (0..YSTART.2) 

CALL  L6AXS  (0..YSTART.IYAXlS.25.YAXIS.90..YMlN.DeLTAY) 
60  TO  170 

160  CALL  L6AXS  (XSTART.O. . IXAXlS.*25.XAXlS.O. tXMIN.OELTAX) 

CALL  PLOT  (XSTART.0..3) 

CALL  PLOT  (0..YSTART.2) 

call  AXIS  (0..YSTART.IYAXIS.25.YAXIS.90..YMIN.DELTAY) 
60  TO  170 

165  CALL  L6AXS  ( XSTART. 0 . . I XAX I S. -25 . XAXI S . 0 . « XMIN » OELTAX ) 

CALL  PLOT  (XSTART. 0. .3) 

CALL  PLOT  (0..YSTART.2) 

CALL  L6AXS  (0. .YSTAPT. IYAXIS.25. YAXIS.90, .YMIN,OELTAY) 
170  INDEXT=3 

60  TO  (190.175.180.185).  ITYPE 
175  YMIN=o. 

OELTAYsl, 

60  TO  190 
180  XMIN=0. 

DELTAX=1. 

60  TO  190 
185  XMlNsO. 

OELTAXal. 

YMlNaO, 

deltay»i. 

190  CONTINUE 

IF  (NFnS.LE.O)  60  TO  195 
FNSM1*NFNS-1 

nXlN»XSTART*OELTAX/FNSMl 
0YINsYSTA«T*0ELTAY/FNSm1 
195  DO  200  Jal.Nl 

X6 ( INOEXT) sX ( J) 

6(INDEXT)=SI6N*Y(J) 

INDEXT»INDEXT*1 
200  CONTINUE 

EPSsEPSl*(ABS(XMIN) .ABS (OELTAX) ) 

N6*N1.4 

X6 ( I ) a-FNSMl*OX IN^XMIN-ABS ( XMlN) -APS ( XG ( 3) ) -1 . 
XG(2)sXG(3)-EPS 


XG(Nl43)*XG(Nl<*2)*ePS 

ZZ*YMIN 

IF  (SIGN.LT.O.)  ZZs-YMIN-50,*OELTAY 

G(1)*ZZ 

G<2)»ZZ 

6(N1*3)»ZZ 

G(N6)aZZ 

CALL  PLOT  (XSTAPT*0..-3) 

X (N1*1)»XMIN 

X(N1*2)*0ELTAX 

Y(N1*1)*YMIN 

Y(N1*2)*0ELTAY 

CALL  LINE  (XfY«Nl*l«0«0) 

OXKKsO. 

DYKK*0. 

PELINCsOELTAX/OELTAY 
XG(N6) sSIGN 
60  TO  355 

205  SIGM=XG(NG) 

XG(NG)sX(Nn 
IF  (NFMS)  225*215»2iO 
210  DXKK=DXXK*DXIN 

DYKKiOYKK^OYlN 
215  DO  220  J*l»Nl 

Y ( J) aSIGN*(Y ( J) ♦DYKK) 

X(J)aX(J)-OXKK 
220  CONTINUE 

225  CALL  LOOKUP  ( X ( 1 ) .XG ( 1 ) . J J) 

IF  ( JJ.GE.MAXDIM)  GO  TO  360 
DO  230  JslfJJ 
XH(J)»XG(J) 

H(J)a6(J) 

230  CONTINUE 

IG*JJ*1 
XH(IG)aX(l) 

H(I6)3F(X(1) *X6(JJ) «G(JJ) «XG(1G) *6(IG) ) 

INDEXGsJJ 

INDEXT=1 

ZlaXd) 

FlaH(I6)-Y(  1) 

IT  = 2 
JJsIG 

IF  (H<16)  .GE.Yd)  ) 60  TO  235 
IF  (JJ.GE.MAXDIM)  60  TO  360 
JJsIG*! 

H(JJ)aYd) 

XH(JJ)sZI*EPS 
235  LASTaO 

X1=Z1 


240 





IF  (XG(16) .LT.X(IT) ) GO  TO  245 
IWHICHsO 
X2sX(IT) 

F2aF(X2.XG(IG-l) .G(IG-l) fXG(I6) * G ( I G) ) -Y < I T ) 

ITaIT4.1 
GO  TO  250 
245  X2sXG(IG) 

IWHICHsl 

F2sG(IG)-F{X2.X(IT-l)  ♦Y(IT-l)  ^XCIT)  fYdT)  ) 

I6aIG*l 

250  IF  (Fl*F2.GT.O.)  GO  TO  260 

OENOMaX2-Xl 

IF  (OENOM.EQ.O.)  DENOMs. 00001 
SLOPE* (F2-F1)/DEN0H 
IGGsIG-1-lWHICH 
ITT=IT-2*IWHICH 

IF  (ABS(SLOPE*RELINC) .GT. l.E-6)  GO  TO  255 

Z2=X2 

60  TO  270 

255  Z2=XI-F1/SL0PE 

60  TO  270 
260  XlsX2 

F1=F2 

IF  (IT.LE.NI)  go  to  240 
265  LASTxl 

Z2aX(Nl) 

CALL  LOOKUP  (Z2 . XG ( IN0EX6) . I6G ) 

IGGsIN0EXG*IG6-1 

ITT=N1-1 

270  ZZ*.99*ZU.01*Z2 

CALL  LOOKUP  { ZZ « X ( INDEXT) .K 1 ) 

CALL  LOOKUP  (ZZtXG ( INDEXG) fKZ) 

K1=K1*INDEXT-1 

K2aK2*lN0EXG-l 

IF  (FfZZ*X{Kl) tYCKl) *X (Kl*l) tYCKl*!) ) .GT.F(ZZtXG(K2) *G(K2) » 
1 XG(K2>I) ,G(K2*l) ) ) GO  TO  300 

IF  (JJ*I6G-INDEX6.GE,MAXDIM)  GO  TO  360 

NGPsITT-IN0EXT*2 

NN2SJJ 

NJJaJJ 

ANUX(NJJ)=XH(NJJ) 

ANUY (NJJ) =H (NJJ) 

IF  (NGR.E0.2)  60  TO  280 
NJlalNOEXT^l 

DO  275  IsNJlflTT 

NJJaNJJ^l 

ANUX(NJJ)aX(I) 

ANUY (NJJ) sY ( I ) 

275  CONTINUE 
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2B0 


285 


2<?0 

295 


300 


305 

310 


315 


NJJsNJJ^l 
ANUX (NJJ)sZ2 

ANUY<NJJ)»F(Z2*X (ITT) .Y ( ITT) tX ( ITT*1 ) * Y ( ITT*1 ) ) 

NNMaNN2^NGP-l 

NL*0 

00  265  IaNN2fNNH 

nl»nl^i 

ANUHX(NL)sANUX (I) 

ANUHY(NL)*ANUY(I) 

CONTINUE 

ANUMX{NL^1)»XMIN 

ANUHY(NL*1)*YMIN 

ANUHY(NL^2)»0ELTAY 

ANUHX(NL*2)*0ELTAX 

IF  (lOASH.EQ.IHY)  CALL  OASHL  ( ANUHXt ANUHY.NL* 1 ) 

IF  (INOEXG.EQ.IGG)  GO  TO  295 

J1*INDEX6^1 

DO  290  I=J1*I6G 

JJaJJ^l 

XH(JJ)sXG(I) 

H( JJ)=G(I) 

CONTINUE 
JJsJJ> I 

XH( JJ)sZ2 

H(JJ)aF(Z2fXG(IGG) «6(I6G) *X6(lG6«l) «G(IGG^1) ) 

INDEXGsIGG 

iNOEXTaITT 

GO  TO  320 

NGRAPHsITT-INDEXT*2 

IF  ( JJ4NGRAPH-1.6T.MAXDIM)  GO  TO  360 
N2*JJ 

IF  (NGRAPH,E0.2)  60  TO  310 
JlalNOEXT*! 

DO  305  laJl.ITT 
JUaJU* 1 
XM(JJ)aX(I) 

H( JJ)aY(I) 

CONTINUE 

XH( JJ)*Z2 

H(JJ)sF(Z2*X(ITT)  fYdTT)  .XCITT*!)  ,Y(ITT*1)  ) 

NM*N2*NGRAPH-1 

L*0 

00  315  IaN2tNM 
L*L*1 

XPTS(L)sXH(I) 

YPTS<L)*H(I) 

CONTINUE 

XPTS(L*1)*XMIN 
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non 


330 

335 


340 


345 

350 

355 


360 


ii».n  i'  ■wf'i^wnwiff'y^ 


XPTS(L*2)*0ELTAX 

YPTS(L^l)»SI6N*YMIN 

YPTS (L*2) »SIGN*DELTAY 

CALL  LINE  (XPTS.YPTS.Ltl*0*0) 

INOEXTaITT 

INOEXGsIGG 

IF  (LAST.EQ.I)  go  to  325 

X1*X2 

Flap? 

ZlaZ2 

IF  (IT.LE.Nl)  GO  TO  240 
GO  TO  265 

IF  (XG(NG) .LE.XG(NG-I) ) NG=NG-1 
IF  (XG(NG) .LE.X(Nl) ) GO  TO  335 
IF  ( JJ43«NG-IGG.GT.HAX0IM)  GO  TO  360 
XH( JJ*l)aXH(JJ) ♦EPS 

JJa J J ^ J 

H ( JJ) aF (X (Nl) »XG( IGG) fG(IG6) ♦XG<IGG^1) ,G( IGG^l) ) 
IGGPlaIGG*! 

DO  330  JalGGPltNG 

JJaJJ^l 

XH(JJ)aX6(J) 

H(JJ)aG(J) 

CONTINUE 

NGaJJ^2 

IF  (NG.GE.MAXDIM)  go  to  360 
DO  340  1*1, JJ 
G(I)aH(I) 

XG{I)=XH(I) 

CONTINUE 

XG(JJ^l)aXG<JJ)^EPS 

G(JJ^l)aYMlN^DYKK 

IF  (SIGN.LT.O.)  G(JJ^l)a-YMIN-50.*0ELTAY^DYKK 
G(NG)sG(JJ^1) 

IF  (NFNS.LT.O)  go  to  350 
DO  345  1*1, Nl 
X(I)*X(I)^OXKK 
Y ( I) *SIGN*Y ( I) -OYKK 
CONTINUE 
XG(NG)*SIGN 
CONTINUE 

CALL  PLOT  (XAXIS»0,0»3) 

CALL  PLOT  (SKIP, 0, 0,-3) 

RETUPN 

MPITE  (6*5020)  MAXOIM 
RETURN 


•♦•♦•format  statements 
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o r>  o r>r>o  o non  o r>  o o r>  r> 


( 


I 5000  FORMAT  (/37H  >>  00  YOU  WANT  ANY  "HI00EN«  3-D  0ATA»/4Xt 

140HPLOTTE0  WITH  A DASHED  LINE  — YES  OR  NO<) 

5010  FORMAT  (/33H  **  INPUT  ERROR  (HIDt):  X{I-1)  s ♦ IPE 1 0 .4 ,5X teHX ( I ) =. 

11X*IPE10.4,/,35H  X(I-l)  MUST  BE  LESS  THAN  X(I») 

5020  FORMAT  (/48H  **  ERROR  IN  DIMENSIONED  ARRAYS  (HIDE):  HAXDIM  »tl4,/, 
135H  INCREASE  MAXOIM  TO  RUN  PROBLEM.) 

5030  FORMAT  <Al) 

C 

t END 

SUBROUTINE  REAOR  (NPASStMASK) 

THIS  SUBROUTINE  HANDLES  THE  INTERACTIVE  I/O  READS  IN  PLOT 

INSTRUCTIONS  AS  WELL  AS  THE  DATA, 

COMMON  /INPUT/  IXAXIS(5) ,IYAXIS<5) »ITITL1 (5) t ITITL2 (5) f ITITL3 (5) t 
1ITITL4(5) ,ITITL5(5) .ITYPE.XAXIS*YAXlS*FACT,ITITLE*Y2fXl*I6PI0f 
2YINC»XINCfNOU«NlN*ILINEStLDEV»NSETSf IHNflHYflHBt IHO 


COMMON  /DATA/  XARRAY ( 900 ) ♦ YARRAY ( 900 ) «HISTOX (900 ) *HISTOY ( 900 > * 
1XPTS(900) ,YPTS(900) «ANUX(900) »ANUY<900) *ANUHX(900) ,ANUHY(900) 


COMMON  /DRAW/  NPTS(20) •IPLOT(20) »IHIST(20) f L INTYP ( 20 ) * INTEO (20 ) * 
ILOGTYP (?0) tSKIP 

...  THE  FOLLOWING  3 CAROS  ARE  FOR  UNIVAC  MACHINES. 

DIMENSION  IA(2) 

IA(l)»6HPAS6fT 
IA(2)>6H  29. 

IF  ( (NPASS) ,GT.l)  60  TO  125 
100  WRITE  (NOU«5060) 

READ  (NINf5030tERR«105)  LOEV 
60  TO  110 

105  WRITE  (NOU»5020) 

GO  TO  100 

...  THE  FOLLOWING  4 CAROS  ARE  FOR  UNiVAC  MACHINES, 

110  IF  (LDEV.EQ.O)  60  TO  115 
GO  TO  120 

115  CALL  EPTRAN  (6«IA) 

L0EV*29 

...  THE  following  CARO  WOULD  BE  USED  FOR  CDC  MACHINES. 

C3  IF  (LOEV.E(J.O)  L0EV*29 
C 

120  CALL  INIT 


46 


I 


125 


130 

135 

1*0 

1*5 


150 

155 

160 


165 

170 

175 

180 

185 

190 

195 

200 

205 

210 


WRITE  (NOU«5070)  NPASS 
WRITE  (NOU«5080) 

READ  (NIN*5050)  MASK 
WRITE  (NOU*5090) 

READ  (NIN.5030»ERRs135)  XAXIS 

60  TO  1*0 

WRITE  (NOU*5020) 

60  TO  130 

IF  (XAXIS. EQ.O)  XAX1SS9 
WRITE  (NOUf5100) 

READ  (NIN«5030*ERRx150)  YAXIS 

60  TO  155 

WRITE  (NOU«5020) 

60  TO  1*5 

IF  (YAXIS. EQ.O)  YAXIS»7 
WRITE  (NOU*5110) 

READ  (NIN.5030»E«Rs165)  FACT 
IF  (FACT.EQ.O.)  FACT»1. 

IP  (FACT.6T.0,)  60  TO  170 
WRITE  (NOU»5120) 

60  TO  160 
WRITE  (NOU»5020) 

60  TO  160 
CALL  FACTR 

IF  (MASK.NE.IMY)  GO  TO  175 
WRITE  (NOUf5010) 

READ  (NIN.50*0)  ITITLl 

60  TO  235 

WRITE  (NOUf5130) 

READ  (NIN.5050)  ITITLE 

IF  (ITITLE. EQ.IHN. OR. ITITLE. EQ.IHti. OR. ITITLE. EQ.IHO)  60  TO  235 
WRITE  (NOU«51*0) 

READ  (NlN,5030.ERRal85)  XI. Y2 
60  TO  190 
WRITE  (NOU.5020) 

60  TO  180 
WRITE  {NOU.5150) 

READ  (NIN. 5030«ERRsl95)  ILINES 
60  TO  200 
WRITE  (NOU.5020) 

60  TO  190 

WRITE  (NOU.5160)  ILINES 
DO  230  lal. ILINES 
60  TO  (205.210.215.220.225),  I 
READ  (NIN.50*0)  ITITLl 
60  TO  230 

READ  (NIN.50*0)  ITITL2 
60  TO  230 

read  (NIN.50*0)  ITITL3 
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220 

60  TO  230 

READ  (NIN«5040) 

ITITL4 

225 

GO  TO  230 

READ  (NIN*5040) 

ITITL5 

230 

CONTINUE 

235 

WRITE  (NOU*5170) 

PEAO  (NIN.50A0)  IXAXIS 
WRITE  {NOU*5180) 

READ  (NIN«50A0)  lYAXiS 


IF  (MASK.EQ.IHY)  GO  TO  255 
WRITE  (N0U»5190) 

READ  <NIN«5050)  I6RID 

IF  (I6RID.EQ.IHN.OR.I6RID.EQ.IHB.OR.IGRIO.EU.IHO)  60  TO  255 
240  WRITE  (NOU*5200) 

READ  (NlN*5030tERRs245)  L6RID 
60  TO  250 

245  WRITE  (NOU»5020) 

60  TO  240 
250  XiNCalO. 

YINC»10. 

IF  (LGRID.LT.O)  XINCsO. 

IF  (LGRID.GT.OI  YiNCsO. 

255  WRITE  (NOU»52lO) 

RFAO  (NlN.5030.ERRs2to0)  ITYPE 
60  TO  265 

260  WRITE  (NOU.5020) 

60  TO  ?55  . 

265  WRITE  (NOU.5220)  NPASS 

READ  «NIN.5030fERR=270)  NSETS 
GO  TO  275 

270  WRITE  (NOU*5020) 

60  TO  265 
275  11=0 

DO  375  Ksl. NSETS 

IF  (MASK.EQ.IHY. AND. K.EQ.l)  WRITE  (NOU»5000) 

IF  (MASK.EQ.I-ir)  GO  TO  360 

280  WRITE  (NOU*5230)  K 

READ  (NIN#5030*ERR=285)  LINTYP(K) 

60  TO  290 

2P5  WRITE  (NOU*5020) 


60 

TO  280 

290 

IF 

(LINTYP(K) .GE 

.0)  GO  TO  295 

IF 

(ITYPE. EO.l) 

IPLOT(K)al 

IF 

(ITYPE. GT.l) 

IPL0T(K)a4 

IF 

(ITYPE. EQ. 2) 

L06TYP(K)=1 

IF 

(ITYPE. EQ. 3) 

L06TYP(K)a-l 

IF 

(ITYPE. EO. 4) 

L06TYP(K)=0 

60 

TO  330 

295 

IF 

(LINTYP(K) .6T 

.0)  GO  TO  315 

o o o o o ooo  r>  o o 


L 


300 

305 

310 

315 

320 

325 


330 

335 

340 

345 

350 

355 

360 


• • • 


WRITE  (NOUt5240)  K 

READ  (NIN*5030.ERR=305)  IPLOT(K) 

60  TO  310 
WRITE  (NOU.5020) 

60  TO  300 

IF  (ITYPE.6T.1)  IPL0T{K)aIPL0T (K) *3 
IF  (ITYPE.EQ.2)  L06TYP(K)=1 
IF  (ITYPE.EQ.3)  L0GTYP(K)=-1 
IF  (ITYPE.E0.4)  LOGTYP(K)aO 
GO  TO  340 

WRITE  (NOU*5250)  K 

READ  (NIN*5030»ERR=320)  IPLOT(K) 

60  TO  325 
WRITE  (NOU*5020) 

60  TO  315 

IF  (IPLOT(K) ,EQ.2)  IPLOT(K)=3 
IF  (ITYPE.6T.1)  IPL0T(K)=IPL0T(K) ^3 
IF  {ITYPE.EQ.2)  LOGTYP(K)sl 
IF  (ITYPE.E0.3)  L06TYP(K)=-1 
IF  (ITYPE.FQ.4)  LOGTYP(K)=0 
WRITE  (NOIJ.5260) 

READ  (NIN»5030»ERR=335)  INTEQ(K) 

GO  TO  340 
WRITE  (N0U»5020) 

GO  TO  330 

IF  (IPLOT(K) .EQ.3.0R.IPL0T(K) .EQ.6)  GO  TO  360 

WRITE  (NOU*5270)  K 

READ  (NIN»5050fERR=350)  IHIST(K) 

60  TO  355 
WRITE  (NOU.5020) 

60  TO  345 

IF  (IHISTCK) .NE.IHY) IHIST(K)sIHN 

IF  (IHIST(K) .EQ.IHY)  WRITE  (NOU»5280) 

WRITE  (NOU.5290)  K 
KFIXsO 

THE  FOLLOWING  LOOP  READS  THE  (X»Y)  DATA, 

DO  365  I=!l»100000 
II=II>1 

UNIVAC  FORMATTED  READ. 

READ  (NIN,5030*ENO=370)  XAHR AY (I  I ) ♦ YARRA Y (I  I ) 

CDC  READ  WITH  2-BRANCH  EOF  CHECK. 

REAO{NIN,1000)XARRAY{II) ♦YARRAY(II) 
IF(E0F(NIN) )54,100 
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o n o 


ClOO  IF{ITYPE.EQ.3.AND,XARRAY(II) ,EQ.O.)KFIX=l 

C 

IF  {ITYPE.EQ.3.AKD.XARRAY(II) ,EQ.O.)  KFIX*! 

IF  {ITYPE.E0.2.AND,YARRAY<II ) .EQ.O.)  KFIXsl 

IF  (ITYPE.EQ,4.AN0,XARRAY(II) .EQ.O.OR.ITYPE.EO.^.AND,YARRAY( 
lin.EQ.O.)  KFIXsl 
365  CONTINUE 

370  1*1-1  , 

11*11-1 
NPTS(K)sI 
375  CONTINUE 

IF  (KFIX.EQ.l)  CALL  FIXUP 
RETURN 

••♦♦•format  statements***** 

5000  FORMAT  (/51H  **  NOTE:  FOP  A 3-0  PLOT  ONLY  A SOLIO-STRA I6HT  LINE./*  i 

141H  IS  available  for  all  data  sets.) 

5010  FORMAT  (/43H  >>  PLOT  TITLE  ( ONE  LINE.  25  CHARACTERS  )<./♦ 

125ri 25  CHAR >) 

5020  format  (/36H  *»  INPUT  ERROR:  RE-tNTER  LAST  LINE.) 

5030  FORMAT  () 

5040  FORMAT  (5A6)  j 

5050  format  (A1)  1 

5060  FORMAT  (/38H  >>  LOGICAL  UNIT  NUMBER  FOR  PLOT  TAPE<»/.  I 

124H  OEFAULT  a UNIT  "29'».)  { 

5070  FORMAT  (/38H  **  THE  FOLLOWING  PERTAIN  TO  GRAPH  NO. .12./) 

5080  FORMAT  {49H  >>  DO  YOU  WANT  A 3-OIMENSIONAL  PLOT  — YES  OR  NO.  | 

U7H  (DEFAULT  = 2-0)0  I 

5090  FORMAT  (/43H  >>  X AXIS  LENGTH  (INCHES)  — DEFAULT  a 9",) 

5100  FORMAT  (/43H  >>  Y AXIS  LENGTH  (INCHES)  — DEFAULT  a 7".) 

5110  FORMAT  (/39H  >>  PLOT  SIZE  SCALING  FACTOR  — OEFAULT. l4H  a FULL  SCA 
ILE.) 

5120  FORMAT  (/43H  **  INPUT  ERROR:  SCALING  FACTOR  MUST  BE  NON. 

110H-NEGATIVE../.23H  RE-ENTER  LAST  LINE.) 

5130  FORMAT  (/5lH  >>  00  YOU  WANT  A TITLE  BLOCK  — YES  OR  NO  (DEFAULT. 

113H  a NO  TITLEX) 

5140  FORMAT  (/51H  >>  GIVE  X.Y  COORDINATE  (INCHES)  FOR  PLACEMENT  OF./. 

14X.37HTHE  UPPER  LEFT  CORNER  OF  TITLE  BLOCK.) 

5150  FORMAT  (/43H  >>  YOU  NOW  HAVE  ROOM  FOR  5 LINES  OF  PRINT.. 

124H  25  CHARACTERS  PER  L INE . . / . 4X . 27HH0w  MANY  LINES  DO  YOU  NEEDO 
5160  FORMAT  (/24H  >>  ENTER  SCRIPT  FOR  THE.I2.9H  LlNE(S):./. 

125H  — 25  CHAP  — >) 

5170  FORMAT  (/38H  >>  X AXIS  LABEL  (25  CHARACTER  LIMITX./. 

125H  — 25  CHAR  — >) 

5180  FORMAT  (/38H  >>  Y AXIS  LABEL  (25  CHARACTER  LIMITX,/. 

125H 25  CHAP  — >) 

5190  FORMAT  (/44H  >>  DO  YOU  WANT  A GRID  — YES  OR  NO  (DEFAULT. 

112H  a NO  GRIDX) 


5200  FORMAT  (/4bH  >>  SELECT  ONE  OF  THE  FOLLOWING  FOR  YOUR  6RID:,/,AX» 
132H  -1*  HORIZONTAL  GRID  LINES  0NLYt/.4Xf 40H  0=  HORIZONTAL  AND  V 

2ERTICAL  GRID  LINES*/f 4X.30H  ♦!=  VERTICAL  GRID  LINES  ONLY) 

5210  FORMAT  (/34H  >>  HOw  ARE  THE  AXES  TO  BE  SCALED< ,/tbX.22Hl  « LINEAR 
IX*  LINEAR  Y*/»bX*19H2  = LINEAR  Xt  LOG  Y»/.6X*19H3  = LOG  X.  LINEAR 
2Y*/,6X,lbH4  a LOG  X*  LOG  Y) 

5220  FORMAT  (/38H  >>  HOW  MANY  SETS  OF  DATA  ON  GRAPH  NO.*I2.1H<) 

5230  FORMAT  (/39H  >>  LINE/SYMPOL  COMBINATION  FOR  SET  NO. * 12, 1H: ♦/fAX. 
147H  Oa  POINTS  CONNECTED  BY  LINE*N0  SY  MBOLS  PRINTED, /♦‘fX, 

248H4N3  POINTS  CONNECTED  BY  LINE,  SYMBOLS  PRINTED  AT,/,4X, 

325H  EVERY  N-TH  DATA  POINT, /,4X,37H-Na  NO  LINE  DRAWN,  SYMBOLS  PR 
4INTED  AT,/,4X,25H  EVERY  N-TH  DATA  POINT) 

5240  FORMAT  (/33H  >>  WHAT  TYPE  OF  LINE  FOR  SET  NO . , 12 , 1H< , / ,6X , 
ll9Hi  a SOLID, STRAIGHT  ,/,6X,20H2  = DASHED, STRAIGHT  ,/,6X, 

216H3  a SOLID, SMOOTH) 

5250  FORMAT  (/33H  >>  WHAT  TYPE  OF  LINE  FOR  SET  NO, , 12, 1H<,/,6X, 

119H1  a SOLID, STRAIGHT  ,/,6X,l7H2  a SOLID, SMOOTH  ) 

5260  FORMAT  (/46H  >>  INTEGER  EQUIVALENT  OF  CALCOMP  SYMBOL  TO  BE,/,4X, 
132HPRINTED  AT  EACH  N-TH  DATA  POINTO 
5270  FORMAT  (/39H  >>  DO  YOU  WANT  A HISTOGRAM  FOR  SET  NO . , 12 , IH, , / , 4X, 
135HYES  OR  NO  (DEFAULT  a NO  HISTOGRAM) <) 

5200  FORMAT  (/50H  **  CAUTION:  FOR  A HISTOGRAM, YOU  NEED  ONLY  SPECIFY,/, 

14X,57H  (1)  A LOWER  AND  UPPER  LIMIT  FOR  THE  FIRST  BIN,  THEN,/, 

24X,54H  (2)  AN  UPPER  LIMIT  ONLY  FOR  EACH  BIN  THEREAFTER,) 

5290  FORMAT  (/35H  >>  ADD  YOUR  (X,Y)  DATA  FOR  SET  NO, , 12, IH; ,/,4X, 

117HEN0  WITH  "PEOF",) 


il 


No.  of 
Copies 


Organization 


DISTRIBUTION  LIST 

No.  of 
Copies 


Organization 


.2  Commander 

Defense  Documentation  Center 
ATTN:  DDC-TCA 
Cameron  Station 
Alexandria,  VA  22314 

3 Director 

Defense  Nuclear  Agency 
ATTN : RATN 
Washington,  DC  20305 

1 Field  Command,  DNA 
Livermore  Division 
ATTN : NIAJ  Brown 
P.  0.  Box  808 
Livermore,  CA  94550 

2 DNA  Information  5 Analysis 

Center 

TEMPO,  General  Electric  Co. 
ATTN:  Mr.  W.  Chaw 
Dr.  Hendrick 
816  State  Street 
Santa  Barbara,  CA  93102 

1 Commander 

US  Army  Materiel  Development 
and  Readiness  Command 
ATTN:  DRCDMA-ST 
5001  Eisenhower  Avenue 
Alexandria,  VA  22333 

1 Commander 

US  Army  Aviation  Research 
and  Development  Command 
ATTN:  DRSAV-E 
12th  and  Spruce  Streets 
St.  Louis,  MO  63166 

1 Director 

US  Army  Air  Mobility  Research 
and  Development  Laboratory 
Ames  Research  Center 
Moffett  Field,  CA  94035 


1 Commander 

US  Army  Electronics  Command 

ATTN:  DRSEL-RD 

Fort  Monmouth,  NJ  07703 

1 Commander 

US  Army  Missile  Research 
and  Development  Command 
ATTN:  DRDMI-R 
Redstone  Arsenal , AL  35809 

1 Commander 

US  Army  Tank  Automotive 

Research  § Development  Cmd 
ATTN:  DRDTA-RWL 
Warren,  MI  48090 

2 Commander 

US  Army  Mobility  Equipment 
Research  5 Development  Cmd 
ATTN:  Tech  Docu  Cen,  Bldg  315 
DRSME-RZT 

Fort  Belvoir,  VA  22060 

1 Commander 

US  Army  Armament  Materiel 
Readiness  Command 
ATTN:  DRSAR-LEP-L,  Tech  Lib 
Rock  Island,  IL  61299 

2 Commander 

US  Army  Armament  Research 
and  Development  Command 
ATTN:  DRDAR-LCN-F,  Bldg.  65 
DRDAR-TSS 
Dover,  NJ  07801 

1 Commander 

US  Army  White  Sands  Missile 
Range 

ATTN:  STEWS -TE-F 
White  Sands,  NM  88002 


aub  ttdMbim 


DISTRIBUTION  LIST 


No.  of 

Copies  Organization 

1 Commander 

US  Army  Armor  and 
Engineering  Board 
ATTN ; STEBB-AD-S 
Fort  Knox,  KY  40121 

2 Commander 

US  Army  Harry  Diamond  Labs 
ATTN:  DRXDO-TI 
DRXDO-NP 

2800  Powder  Mill  Road 
Adelphi,  MD  20783 

1 Director 

US  Army  TRADOC  Systems 
Analysis  Activity 
ATTN:  ATAA-SL,  Tech  Lib 
White  Sands  Missile  Range 
NM  88002 

2 Commander 

US  Army  Nuclear  Agency 
ATTN : MONA-WE 
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